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1. Introduction a 

To describe mathematically the mechanical behaviour of a continuous medium x 
it is first necessary to define tensors specifying the stress, the deformation and (4 
- any required time rates of variation of these two quantities. These tensors must 
be defined in such a manner that the necessary requirements for invariance under c- 

_ rigid body motions are satisfied. This problem has been studied from various 5 
; points of view by OLDROYD (1950), TRUESDELL (1955, a, b, c), RIVLIN & ERICKSEN ii 
(1955), NoLt (1955), GREEN & RIVLIN (1957) and others. The mechanical prop- * 
erties of the medium are then described by means of functional relations between 
the mechanical and kinematic tensors thus defined. om 
An important class of problems arises when one of the mechanical tensors 

is connected with the remainder by means of a polynomial relationship. To be “ 


explicit, we may suppose the stress tensor to be expressed as a polynomial in 
the elements of the ~+1 tensors defining the deformation and its first m time 
rates of variation. The problem then reduces to the determination of a matrix 
polynomial suitably restricted to describe any symmetry properties of the 
material. When the body is isotropic the stress matrix may be expressed as the 
sum of a finite number of terms, each term being a matrix product formed from 
the kinematic matrices, multiplied by a scalar invariant function of these ma- 
trices. Each of these invariant functions may in turn be expressed as a poly- 
nomial in a finite number of invariants. This reduction has been carried out 
by Rivitn (1955) for a polynomial in two matrices and by SPENCER & RIVLIN 
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(1959a, b) for the general case of m matrices. The writer (1958) has considered © 


the corresponding problem for transversely isotropic and orthotropic bodies when 
two kinematic matrices are involved. 

In the present paper, this work for aeolotropic materials is extended to the 
general case involving m matrices, attention being concentrated mainly on the 
transversely isotropic case which is treated from a slightly different viewpoint 
from that previously adopted. A preliminary reduction is achieved using results 
for the invariants of linear and quadratic forms under transformations of the 
orthogonal group (WEITZENBOCK 1923, TURNBULL 1928). The results are sim- 
plified by using the matrix relations derived by SPENCER & RIvLIN and by 
further applications of the Cayley Hamilton theorem. For transversely isotropic 
bodies an alternative formulation is possible based on the decomposition of the 
3 x3 symmetric kinematic matrices into a system of symmetrical and unsym- 
metrical 2x2 matrices. This formulation is given in § 6 and relations between 
the two alternative systems of invariants for two matrices are derived in § 7. 
These relations are readily obtained by using the symbolic method of classical 
invariant theory, and by an extension of this method, given in § 8, the total 
number of coefficients and invariants required for the description of transverse 
isotropy may be considerably reduced. The stress strain relations for an elastic 
body resulting from the two formulations are given in §9. The generalization 
to nm +1 matrices for orthotropic bodies is indicated in § 10. 


2. Notation and formulae 
We consider a body of homogeneous material each particle of which is at 
rest at time ¢=G and which subsequently undergoes a continuously varying 
deformation. A system of curvilinear coordinates # is associated with elements 
of the material and moves with the body as it is deformed. The covariant and 
contravariant metric tensors for this system at the initial time t=O are 
¥i;=7:;(0",0) and y'’=y"!(9", 0) respectively, the corresponding quantities at 
the current time ¢ being Jj, =J},(0’, t) and I'’=I"’ (8, t). In terms of these 

quantities the covariant strain tensor »,; is defined by 


1 
ij = 2 (Lig — Vig), (2.1) 
and by differentiation we may derive successive covariant rate of strain tensors 
af), al?) ay, aft. Thus 
Clair URN ee ot oe EF) Ti 
rag iat pacietae (Pier D eee) | (2:2) 


where D/Di denotes differentiation with respect to ¢ holding the convected 
coordinates # constant. Expressions for these quantities referred to frames of 
reference fixed in space have been given by OL_pRoyp (1950) and others: 
Riviin & ERICKSEN (1955) have derived equivalent expressions for the descrip- 
tion of the motion of a deforming material, again referred to a fixed frame of 
reference, their starting point being the displacement, velocity and successive 
acceleration gradients. 

In the present work it is convenient to use the metric tensors vit, y;; to raise 
and lower affixes. We therefore write 


i= n;, nisy'*n, (2.3) 
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with corresponding expressions for the remaining kinematic tensors derived from 
(2.2). When the curvilinear system # is such that it coincides with a rectangular 
Cartesian system x’ (= x,) in the undeformed body at time ¢=0 we denote the 
quantities 7; ;, a”) by @;,, af respectively. With this choice of coordinates the 
metric tensors y,;;,y'’ reduce to Kronecker deltas 6;;, 6‘ and the covariant, 
mixed and contravariant components derived by the process (2.3) are equal. 
From these components we form the symmetric kinematic matrices 


E=|le;; 


|, A, =||a?)||, (Ae eae (2.4) 


The contravariant stress tensor referred to the convected system 9 is denoted 
by 7’, or by ¢’’ if #&=x, when t=0. A homogeneous, uniformly aeolotropic 
material may then be defined by the tensor relation 


T= PC, 6 Mpg Am) =P"), (2.5) 
where the functions /"’ are polynomials in the arguments indicated. Any coeffi- 
cients occurring in these polynomials are material constants which are independent 


of position throughout the material and do not depend upon the deformation. 
More generally, in place of (2.5) we may write 


ee de fA, (2.6) 


where ie are symmetric polynomial tensor functions of the kinematic tensors 
Gray al) (k =1, 2,...,) while the quantities gw, are invariant functions of these 
tensors which are not necessarily polynomials but which are without singularities. 


4 


By a suitable choice for the functions 9,, ioe (2.6) may be made to yield the 
stress deformation relation for aeolotropic bodies derived from a different point 
of view by GREEN & RIVLIN (1957). We shall here be concerned mainly with the 
form of the functions ie and shall therefore restrict attention to the simpler 
expression (2.5) from which, by a tensor transformation we obtain 


5 08 06 
oi ree 1S 
4 Ost! Te ? 


ie 
=e Ai. fs = Gigs, ( 7) 
where 
* ae a0 ij 1 ij ij 
A= Oee! Gee C= 2 (Ay, sis Ay) c (2.8) 


Let Q be an arbitrary symmetric 3 x3 matrix with elements ¢;,; referred to 
the coordinate system x;. If we form the quantity 


Gy, Cry, 07) oped es (2.9) 


then F is an invariant polynomial function of the arguments q;;, é,;, a’) which 
is linear and homogeneous in the elements g;;._ Also, is symmetrical in the 
components g;; and in the components of the kinematic tensors, so that inter- 
change of any pair of suffixes of g;; leaves the form of F unchanged. Further- 
more, from (2.9) we have 


ti} a 


1 ( oF OF et 


DAO 
2 \ Oi; O”Fji ( ) 
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Symmetry properties of the material are investigated by considering the 
behaviour of the relations (2.5) under rotations and reflections of the x,-axes 
in the undeformed body. Under all transformations of the orthogonal group 


sa Naess %,=M,;x;,  Mi,Mj,=“,,M,;=5;;, (2.11) 
we have, for each of the tensors 9; ,, ¢;, ay? ,t , t'? the result 
463 = Mi Ms Mrs (X25 = Gar Cp @ ij Dl Gime d, Deane (2.12) 


If the relations (2.5) are form invariant under all transformations of a given sub- 
group of the orthogonal group (2.11) then the function F is also form invariant 
under this subgroup of transformations. 


3. Some invariance relations 
For subsequent developments we require results in classical invariant theory 
which have been given by WEITZENBOCK (1923), TURNBULL (1928) and others. 
Tet aS 
h=4j Vi; fe = O5G Vij 0 000s Im = Mag Vi; (2,7 =4, 2, 3), (3.4) 
denote a set of ternary quadratic forms in the variables y; where, here and sub- 
sequently, the summation convention is used for repeated indices. The forms 
hi» fo -+- fy Taay be represented in symbolic notation by the expressions 


h=(4%) (@¥;) = (a; y;) (a ;¥;) = (a; ¥;) (a;’ ¥;) oe 
fa = (8: i) (8; 4) a ¥) (5; ”) = (o vy) (8 mar 


(3.2) 
hn = (™; y;) (m; Vi) = On yi) (m; . = bn! ys) ont Yj) So, 
where 
4; = 4, =0,4—=4,4,—4;4 =, (3.3) 
and corresponding formulae hold for 0;;,..., m,;;.. Any product expressed in 


symbolic notation then only has any significance if the symbols a;, a;, aj, ..., 
m,,m;,m;... occurring in a particular representation of a given form occur in 
pairs. The use of several symbols a;, a;, a... to represent the same form enables 
any product of the symbols to be expressed uniquely in terms of the coefficients 


a;;,6;;,...,m,;;. For example 


a, 4,4; a; = 4; 5,4; ;, A; 4; 4; A; = A; ;4;;. 


We now consider the behaviour of the forms (3.1) under the group of trans- 
formations 


w=, Vo = MA oo, (3.4) 
where the quantities M,, satisfy the relations 
MoM, = IE Ma Ory glee ts (3.5) 


and, here and subsequently, greek suffixes are restricted to take the values Pape 
while latin indices cover the full range 1, 2,3. The transformations (3.4) form a 
subgroup of the proper orthogonal group 


¥; = M59; (M,; M;, =M;;M,;=6 


M,\= a, 4-4). 


jk? 
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Under the group of transformations (3.4) the forms (3.1) become 
h=4;; 9:9; xb Nig tees bn =M;5 Vi 0; 


where the coefficients @;; are given by 


Big = M,,, Mor Any» ay, =M,,%y)> a,=4)1, (3.6) 
with similar expressions for bij. ...,M,;;. For the symbols a,, a;, a;,..., ™;, 
m;,m; ... we have the transformations 

a; = ay, a, =M,,4,, (3a/) 
corresponding to (3.4). 
Combinations of the coefficients a;;, b;;...m,;; Which are invariant under the 


transformations (3.4) may be derived by adjoining to the quadratic forms (3.1) 
the linear form 


“Ope, (3.8) 
where 
(4, 4, 43) = (1, 0, 0). (3-9) 
It then follows (see, for example, WEITZENBOCK 1923, chapter X) that any 
polynomial invariant formed from the coefficients a;;,b;;,...,m,;; may be 
represented symbolically as a sum of products of the quantities 
C1 ea es 
m Ne Ms\=(Enc), §9; = (En), (3.10) 
G , Cs 
in which £;,7,;,¢; may represent any of the symbols a,, aj, a;,...,m;,™;, Mj; ... 


or the set of coefficients /; of the linear form (3.8). Alternatively, in the symbolic 
expression of any polynomial invariant, the system (3.10) may be replaced by 


(En/), Ey Na = (E|), €,l,= (él) =&, (3.11) 


if we remember that 
(7) = (E|m) + (E4) (2), 
(&C) = (El) (06) — (E02) (2) + (nS) (26). 
If we consider polynomials in the coefficients of (3.1) which are invariant not 


only under transformations of the subgroup (3.4), but also with respect to the 
reflection 


(3.12) 


(V1, Vo Vs) = (V1, Yo, — Vs), (3.43) 


the determinants (£7 ¢), (7 1) may be omitted from the symbolic representations 
(3.10), (3.11). For since such determinants change sign under the transformation 
(3.13), an even number of them must occur in each symbolic product, and the 
product of each pair of these determinants can be expressed as a polynomial in 
the scalar products (£7) or (¢|7) and (J €). It follows that any polynomial in the 
coefficients a;;, b;;,..-,™;; which is invariant under the group of transformations 
defined by (3.4), (3.13) can be expressed symbolically as a polynomial in the 


Sepa ee ee) a (3.14) 
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or as a polynomial in the products 


(Ely) =ata, (6) =4:5;= 41, (3.15) 


ut 


in which £,7 may represent any of the symbols q;, dj, aj, ..., ™j, M;,M; ,.... 
Evidently, if the polynomial basis is abandoned, the systems (3.14), (3.15) are 
sufficient to represent the invariants of the transformations (3.4) alone. 

When the symbolic invariants are expressed in terms of the actual coefficients 
of the quadratic forms (3.1) it follows from (3.14) that any polynomial invariant 
formed from these coefficients can be expressed as a polynomial in the quantities 


As Wie Sim Panes (i) 
Ag 5 bag - Ets Ponts (ii) 


where A; ;,4;;,--+»&;;, ¥;; may represent any of the sets of coefficients a, ,, Oj 5.--M3 
no restriction being placed, for the present, upon the number of times in which 
each of these coefficients may occur in the products (3.16). Alternatively, from 
(3.15) we may infer that any polynomial invariant may be expressed as a poly- 
nomial in the quantities 


(3.16) 


Aug Haan Baa Ves (i) 


& (3-17) 
Wi oo oi: (11) 


4. Transversely isotropic materials 


A material is defined as transversely isotropic with respect to the x-direction, 
if the relation (2.5), or the function F given by (2.9), remains form invariant 
under the group of transformations of the coordinates x; in the undeformed body 
defined by 


t= %, | %y=M,,x 


Q0~oO? 


(4.1) 


(%1,%y,%g) = (%, Xe, +3), 


where the quantities M,, satisfy the orthogonality conditions (3.5). Since (4.1) 
combines the transformations (3.4), (3.13), and the tensors q;;, e;;, a, t? of §2 
behave under these transformations in the same way as the tensors a;;, b,;, ...,;; 
of § 3, we may infer that / may be expressed as a polynomial in the products 
(3.16) or (3.17) where 4;;, 4;;,---,&;;,%;; NOW represent elements of the matrices 
Q, E, A, and the resulting products are linear in g;;._ We consider first the results 
which follow from the forms (3.16) and defer consideration of the alternative 
formulation (3.17) until § 6. 

Remembering the definition (2.9) for F and the symmetry property (3.16) 
we may ess 


M 
F=45 (02 97 On + OR: G;, 70 Qi] +a 24 [Regd hes. (4.2) 
In this expression, Ps, Ope OR. are elements of matrix products of the form 


P=||P,;|j, Es j= O54 Cay. Be Pe (4.3) 


where b,;, ¢;;,--.,4,;,/;; may represent elements of any of the kinematic matrices 
E, A, or may be replaced by the Kronecker delta 6;;; special forms of these 


Curvilinear Aeolotropy 199 


coefficients are obtained by replacing one or more of ” YO or RS by 0;;- 
The functions ®,, Y% are polynomials in the invariants 
[,= = dF c : Amn tna 
ae (4.4) 


de 7 by; Cikes: Ching wi drias 
The quantities g;; are absent from the products One Oy. MR 
variants I,, j,. 


From (2.10) and (4.2) we obtain 


- N M 
= 8 TPR On + OR OOnG +E LR +R IB, (45) 
= ei" 


,; and the in- 


and by a tensor transformation using (2.7) and (2.8) 
N 


mM 
i = OB Ci 0, , + 2 Chis Rg B. (4.6) 

The invariants J, and the coefficients of Y% in (4.6) are the forms appropriate 
to isotropic pee and may be expressed in more familiar forms in terms of the 
tensors 7; ;, as” and the metric tensors J},, I, These isotropic terms have been 
studied by SPENCER & RIVLIN (1959a, b) who have derived integrity bases for 
the invariants and have effected reductions of the symmetric matrix polynomials 
by using the Cayley Hamilton theorem and Peano’s theorem. We are here con- 
cerned with the additional terms and invariants which arise from the assumption 
of transverse isotropy. In the expression for t'’ the quantities Ci, obtained by 
giving v and s the values 1, 2, 3 form a symmetric 3 X3 matrix which exhibits 
the same behaviour under orthogonal transformations of the x; coordinates as 
the matrices H, A,. Denoting this matrix* by C and defining the matrices 


(7) 
Q by ()p — tae A 99S] (4.7) 
we may write AO Oe (PCO), (4.8) 


where the symbol { ],; shane that only the leading element of the matrix 
product is required. The expression (4.8) is therefore of the same form as those 
for the transversely isotropic invariants j, but involves, in addition to the 
kinematic matrices E, A, the transformation matrix C. 

The reduction to simpler forms, both of the transversely isotropic coefficients 
and the invariants J, thus involves a consideration of the leading term only of a 
matrix polynomial of symmetric matrices. Since from (4.3) 


Py a by 5 C3; es dinn Tat = hi q;; + Onn Dnt , (4.9) 
we may equally well regard P,, as the leading term of the symmetric product 
ee a > [O53 Cij aiake dinn Jas = bs Cri Ace Ginn Pale (4.10) 


x Although there are in fact six matrices C obtained from the six independent 
choices for 7,7, symmetry properties are examined by considering transformations of 
the x; coordinates (represented by the suffixes 7, s) in the undeformed material alone. 
The indices i,j preserve the tensor character of t'’ under transformations of the o 
coordinates. We shall therefore refer to C as if only one matrix were involved. As 
an alternative procedure, we could consider the function F and derive the stresses 
by differentiation after the reduction of F has been carried out. The matrix Q would 
then replace C in the analysis. 
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5. Some reductions of symmetric matrix products 
By applications of the Cayley Hamilton theorem and Peano’s theorem, 
SPENCER & RIVLIN (1959a) have shown that any symmetric isotropic matrix 
polynomial in 7 symmetric 3 X 3 matrices a, (P=1, 2, .-., v) may be expressed 
as a linear combination of the matrix products 


Tl a, - @, (I) 
a,ajtaja;, a,a7+a7a;, 4G; a? + a; ai, (II) 
a,4;4,+4,4;4,, a?a,a,+a,4,a;, a;aja,+a,a;a;, (IL) 
aaza,ta,ara?, (i) a,aaja,+a,a;a,a,, (ii) av) 
a, 4; 4,4; + a; a@,a;a;, (iti) 

a; A; A, A, + A, &, a; a;, (V) 

a? a; @,a,+ a,4,4;,07, 4,47 a,a,+ a), 45 4;, (VI) 6-1) 
a? @ a, a, + 4,0,07 07, a, a7 a, a, + a, a; a, (VII) 
a,4;4,0,0;-+ a5 a,4,0;4;, 4; 4; , 0,4) + a) Ga, a; G, (I) 
a; @, A; A, a; + a;a,4,4,4;, 4,4; a, a} a, + a, a? a, a; 4,, 

@; A; A, A, A, + A, A, , A; A; , (IX) 
a; 4; , 4,05, + 4,0, 4,0;,4;, 4; 4, 4, a, a; + GG), @, @, a;, (x) 

a, 4,4, 4; 4;+ a; 4; 45, a, A, 

a; 4; A, A, A, A, + A, Ay, A, A, A; Q;. (XI) 


In the table (5.1) all possible selections of the suffixes 7,7, k, J, m, m (all different) 
are chosen from the numbers 1, 2, ...,7; the coefficients in the resulting matrix, 
polynomial are polynomial scalar invariants under the full orthogonal group 
that is, polynomials in the first kind of invariant J, in (4.4). 

For transversely isotropic bodies, the additional invariants J, are all included 
in the set obtained by taking the leading terms of the products (I) to (XI) and 


substituting the matrices E, A, (p=1, 2,..., ) for @,,a@,,...,@,,,. The coeffi- 
cients (4.8) are obtained similarly from the corresponding set for +2 matrices. 
In this case we substitute C,H, A, (6=1,2,...,m) for a,,@5, ..2, 4,4, and 


choose only those invariants in which C (=a,) occurs linearly. 

By further analysis based on the Cayley Hamilton theorem the writer (1958) 
has shown that the products (ii) and (iii) of set (IV) can be expressed in terms 
of each other, together with products of lower degree in the matrices a,; either 
of the types (ii) or (ili) may therefore be omitted from the scheme (5.1) without 
loss of generality. Furthermore, from a systematic examination of the products 
(VII), (VIIT), (X) and (XI) carried out by SPENCER & RIVLIN* we may infer that 
these may also be expressed in terms of products of lower degree and may thus 
be omitted from (5.1). In the resulting scheme, the symmetry of the system is 


* Private communication. I am indebted to Professor A. E. GREEN for drawing 
my attention to this work which, I understand, is to be published shortly. 
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preserved at the expense of still further redundancies in the coefficients. For 
example, from any four matrices a;, a;, @,, @, we may derive twelve symmetrical 
products of type (V) by taking all possible permutations of the letters i, 4, k, l. 
Not all of these products are linearly independent, however, for relations between 
them follow as a consequence of the Cayley Hamilton theorem. We may observe 
that any reduction in the number of coefficients (5.1) leads to a corresponding 
reduction in the number of invariants J, and J. 


To achieve these further reductions we require the Cayley Hamilton theorem 


a> — a@ tra+ 3a|(tra)?—tra?] —41[(tra)? — 3 tratra?+ 2tra®] =0 


(tra = tracea), Oe 
for a 3 X3 matrix a, and the generalization 
abc+acb+bac+bea+cabtcba 
=(be-+cb)tra+ (ac+ ca) trb+ (ab+ba)tre+ 
+a(trbe —trbtre) + 6(trea — tretra) + 5.3) 


+ e(trab —tratrb) + 
+I(tratrbtre —tratrbe—trbtrea — 
—tretrab+trabc-+treba), 


for three matrices a, b, e obtained by RIVLIN (1955) by substituting a+b and 
a+6-e for a in (5.2). 

In each of the equations (5.2), (5.3) the left hand member is a matrix poly- 
nomial of the third degree in a, b, e while all of the terms on the right hand side 
are of lower degree, the coefficients of these lower degree terms being scalar 
polynomial invariants. To abbreviate the subsequent analysis we shall use the 
following convention. If a matrix polynomial P which is homogeneous and of 


the 7 degree in matrices a, b, c,... can be expressed as a polynomial in which 
all of the terms are of degree <7—1 we shall describe P as reducible and write 
P=0. 


With this convention, (5.2) and (5.3) may be re-written 
a= 0, (52°) 


abc+acb+ bac+bea+cab+cba=o. (5.3°) 


If P and Q are two matrix polynomials, homogeneous and of degree 7, for which 


ky P +kh,Q =0, (5.4) 
k, and k, being numerical coefficients, then obviously we may write 

PPh, O: (5.5) 
In general, if PB, P,,...,P, are a system of m matrix polynomials which are 


homogeneous and of degree 7, and 


a,,P+o.;R+---+4,,BR =0 C= Apes > <1) ; 
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are m linear relations between them with constant coefficients «7, then m of 
the matrices P can be expressed linearly in terms of the remainder together 
with terms of degree <r~—1 by means of the ordinary processes for solving 
simultaneous equations, provided the coefficient matrix ||«,,|| has rank mi. 


In dealing with symmetric products we use the notation 
(abc...def)=abcec...def+fed...cba, (5.6) 
for symmetric matrix polynomials and write 
[abe...def], (5.7) 
for the sum of all terms formed from the product abc...def by a full per- 
mutation of the letters a, b, c,...,d, e, f. With this notation, (5.3) or (5.3’) may 
be written as (abc) + (acb) + (bac) =0, 
Se as [abc] =0. (5.8) 


We define the symmetric multiplication of the symmetric product abc...def 
by the product fu...vw as 


tu...vwiabc...def)+(abe...def)wv...ut 


5.9 
=(tu...vwabc...def)+(abc...defwv...uf), oo 


and observe that a symmetric polynomial results from this operation. We define 
the symmetric substitution of d by tu...vw in the symmetric product 
(ab...cde...fh) by 
(ab...c(tu...vw)e...fh) 
=@0...ctu...0we..fhkthf...etu. twee... 0a-- 54 
10 
+ab...cwv...ute...fh+hf...ewv...utec...ba ( 
=(ab...ctu...vwe...fh)+(ab...cwv...ute...fh), 


the result again being a symmetric polynomial. 


From any four different matrices a, b, c,d we may form the twelve sym- 


metric products 
(abcd), (abdc), (acbd), (acdb), (adbe), (bacd), (i) 
(adcb), (bade), (bead), (bdac), (cabd), (cbad), (ii) (5.11) 


of type (V). Between these, we may obtain by the symmetric substitution of a 
by ad in (5.3’) or (5.8) the relation 


(ad be) + (adcb) + (bade) + (cad b) + (bcad) + (ebad) =0, 


gh [(ad) be] =0. (5.12) 
By similar substitutions of this type, as indicated, we may derive the five further 
relations [(ab) ed] =0, [(ac) bd] =o, 

[((be)ad[=0, [(bd)ac]=0, (5.13) 


[((cd)ab] =0. 
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These six relations are linearly independent in the products (5.11) and yield the 
solutions 


(adcb) = — {((abed) + (abde) + (acbd) + (acdb) + (adbe)}, 
(badc) = (acbd) + (ad be) — (bacd), 

(bcad) = (abdc) + (ad be) — (bacd), 

(bdac) = (bacd) + (abed) — (adbe), Cae 
(cabd) = — {(abed) + (abdc) + (acbd) + (acdb) + (bacad)}, 
(cbad) = (bacd) + (acd b) — (ad be). 


Each symmetric product of the set (ii) in (5.11) can therefore be expressed as 
a linear combination of the products (i), together with terms of lower degree. 
We observe that in (5.14) the formula for (ade b) or for (ca bd) may be 
replaced by one of the relations ; 


(adcb) — (cabd) = (bacd) — (ad be) 
= (acbd) — (bade) 
= (abdec) — (bcad) (5.15) 
= (bdac) — (abcd) 
= (cbad) — (acdb), 


connecting (adc b) and (ecabd). No further reductions result from the sym- 
metric multiplication of (5.8) by the remaining matrix d, for it can be shown 
without difficulty that 


3 (d[abe}) =3d[abe] +3[abel|d 
= 2[(be) ad] + 2[(ae) bd] + 2[(ab) ed] — (5.16) 
— [(ad) be] — | (bd) ac] — [(ed) ab], 


with corresponding formulae for (b[acd]), (e[abd]) and (a[bed}). For any 
given choice of the suffixes 7,7, k,/ from the numbers 1 to 7 we may therefore 
obtain the set of products 

(4;4;4,4)),  (@;4;4,4,), (4, 4,4;q)), 


(5.17) 


(@;@,4,4;), (@;@4;4,), (4,4; 4, @)), 


of type (V). An equivalent result may be obtained by differentiation from the 
integrity basis for invariants involving five matrices derived by SPENCER & 
RIVLIN (1959b). 

Analogous reductions may be carried out for the remaining higher degree 
products of types (VI) and (IX), and a large number of redundant terms eliminated 
in this way. It is, in fact, just possible that products of type (IX) may prove 
to be completely reducible. For if we denote the system (5.14) by {a bed} and 
the sets of six equations derived by the symmetric substitutions of @ by ae, 


= 
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(ae) 

(ae) bed}, fa(be) cd}, fab(ce)d}, {abc(de)}, 
(ad) bce}, {a(bd) ce}, {ab(ed) e, 

(ac) bde}, {a(be) de}, {(ab) cde}, 


yield exactly sixty relations between the sixty symmetric products derived from 
(abede) by permutations of the letters a,b,c,d,e. If these relations are 
linearly independent, each of the products of type (IX) can be expressed entirely 
in terms of products of lower degree. We shall not, however, examine these 
further reductions here, since for the transversely isotropic case all of the in- 
variants I,, J, formed from matrix products of higher degree than the fourth 
by taking the traces and the leading terms respectively can be reduced completely 
by the method of § 8. 


6. Transversely isotropic bodies: an alternative formulation 


An alternative description of the invariant system for transversely isotropic 
bodies may be obtained from (3.17). We again consider the system of 7 sym- 
metric 3 X3 matrices @,, @,,...,a, but in this section we shall denote by b, the 
2x2 matrix derived from @, by omitting the first row and column, and from the 
remaining elements of the matrices a, construct an additional system of 2x2 


matrices. Denoting any two different matrices a;, a, of the set a, (pb =1, 2,...,7) 
for which 7<k, by a, b respectively, so that 

Bi feo a=lla;;||, b= (2, (6.1) 
b; = aval ) C; == || 21 ag1| , Crp poo || 21.4 bex|| (7 = k, Dx = die, a, B = 2, 3) i (6.2) 


This gives a system of 27 symmetric 2 x2 matrices b;, ec; and $7(y —1) unsym- 
metric 2x2 matrices ¢;,. From (3.17) it follows that any of the isotropic or 
transversely isotropic invariants for the system of 3x3 matrices a, may be 
expressed in terms of the isotropic invariants for the system of 2x2 matrices 
b;, €;, ¢;, together with the leading terms (e.g. a,,) of each of the matrices a,. 
Furthermore, each of these invariants of the system (6.2) can contain at most 
one of the matrices ¢€;, €;, linearly. 

To reduce the system of invariants derived from (6.2) we make use of the 
following results derived by RIvLIN (1955) by applications of the Cayley Hamilton 
theorem: 

(7) Any matrix polynomial P in symmetric 2 x 2 matrices B, (p =1, 2, ..., ”) 
can be expressed in the form 


| Pol + LOBit 2 2 Pix B; B, (1 =), (6.3) 
= i= =I 
where g, y; and y;, are polynomials in the invariants formed from the ma- 
trices B, and I is the unit matrix. 


(72) Any polynomial invariant formed from the elements of the system of 
2x2 matrices B, can be expressed as a polynomial in the quantities 


iB; itr Be Bx. 
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From (3.17) we observe that any polynomial invariant of the matrices (6.2) 
is expressible as a polynomial in the traces of these matrices. It follows im- 
mediately from (7) that all invariants formed from the symmetric matrices 
b;, ©; are expressible as polynomials in the quantities 


tO 4tt 0D, tt Cet OG, (6.4) 


if we remember that all invariants must be linear in e;. Again, any invariant 
containing elements of ¢;, must be expressible in the form tr e,;,P where P is 
a matrix polynomial containing only the matrices b,. With this result, the ex- 
pansion (6.3) implies that any polynomial invariant containing the matrices ¢;, 
is expressible as a polynomial in 


tre;,,  tre,;,b,, tre,;,b,b,, trb,;, tr b;b; (2 = m). (6.5) 


Remembering (6.4) and (6.5) we deduce that the invariants for the trans- 
versely isotropic case can be expressed as a polynomial in the quantities 
A1) Deo» Ae p Op x» Ayo 1s Aya xp Opi» Aya Cap 4py yr, (6.6) 
where 41, 0x, Cxp»4,g represent elements of any four of the matrices a; subject 
only to the condition that in the last of (6.6) c,g, d,, are elements of different 
matrices and likewise the matrices a;,;, b;; are not identical. If c,,=d,, the last 
of the invariants (6.6) can be expressed as a polynomial in the simpler invariants 
by using the Cayley Hamilton theorem for a 22 matrix. The system (6.6) 
represents a different expression of the total set J,, J, of the invariants (4.4) and 
not merely of the invariants /j,. 


7. Invariant relations for two matrices 


The connexion between the formulation of §6 and that of the preceding 
sections may be illustrated by examining the systems of invariants formed from 
the two symmetric 3 x3 matrices 


a=|la;; 


b= [ill (7-1 
Corresponding to the theory of §§3 to 5 we have the system 
lo=@,. +1,=-0,,¢ 


aq"G 4? 
Ii=b I, =b,,0 


ig -gu 


T3 = 4; 54; 4,% 5» y=, 15 =; 41, 

at? Ts = 0; 50; On: T= be Ts = b,;6;1, (7.2) 
= a; 5955, Jg= A; 4; Ons, I3= A; 5054 Ons» I= i; jp On O15, 
I3=%Gi0, Se= Ay 4; 50,1, = Ay 555041, Ig = Ay 455954 nr: 

Of these, the invariants J,, I,, 13, 11, 12, 13,1, Jo, Js, Ja are those appropriate 

to isotropic bodies and have been given by RIvrINn (1955). The remainder apply 

only to transversely isotropic materials. From (6.6) we obtain the alternative 

set of invariants 

Ky = a4, Ky = Ay 48 q; K3 = 41, Ky =% e415 Ks = %44.%8 41, 

VS ies Ne Dap Op a» Keg Og Ce Ob Fg = Og Oye Onis is 

Ly = 4p O60; Ly =o 0a, Lg = 4x5 Op1 Ly =e bu9 O61; 


Ls= yg xg %p1» Le= 015 4x8 %81, Ly = M4 44x 5 Op 0,1, DPg= 4 404848) 0y1- 
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To express the system (7.2) in terms of the invariants (7.3) we note that for 
any 2X2 matrix ¢ we have the Cayley Hamilton theorem 


e=ctre+ iI[tre— (tre), (7.4) 
and by post multiplying by any other 2 x2 matrix d and taking the trace of the 
resulting equation we deduce that 

tre’d =tretred + dtrd|tre?— (tre)?]. (7.5) 


Again, if we multiply corresponding sides of (7.4) by the analogous relation for d 
and take the trace of the resulting expression we have 


tre’d? = i tre® trd*+ tredtretrd — (tre)? (trd)?. (7.6) 


Relations between the invariants (7.2) and (7.3) are now readily obtained by 
translating these into the symbolic notation of § 3, employing the first of (3.12) 
and using (7.4) to (7.6) where necessary to simplify the expressions obtained after 
regrouping the symbols 4a;, aj, a;,...,0;,0;,6;,... into the elements a,,, 0;;. 


BAT IE 


This process yields, in the case of J,, for example 
Jy = 4; 5% p Op = (4) (@’B) (64) 
=[(ala’) + (La) (La’)] [(a’ |b) + (La’) (20) | [(b| a) + (22) (Za) 
= Deg Ag y Og + Dy Ong Agy + 2.) 4 Ay g Ops + 
$2414 Part O11 Ma art 41 br 
The first of these terms can be simplified by making use of (7.5) with a, b for 


c, d respectively; the remainder can be recognised as invariants of the set (7.3). 
By means of this procedure we obtain the relations 


hie 
I, = K,+2K,+ Kz, 
I, = 3K, K,— 3 Ki+3K5+ 3KyKy+ Kz, 
[p= Ks, 
T,= K,+ K, 
ti Lye 2L,+ K, Ks, 
Jg= 215+ L5+ K,L,+ 2K,L,+ } K,(K,— Kj) + K3(K,+ K3), 
Jg=2L,+L1,+KiL,4 2K3L,+ ¢ K, (Ky — K,*) + K,(Ka + Kg’), (7.7) 
Ju = 21,4 134 2(KgL,+ K3L,) + 2K,K3L,+ Kiki L,+ 
+ K,L,+ K,L;+ (K,+ K2) (Ky + Ky?) PEG KL RK yee 
+ 2 Ky (Ke — Ky") + 3 Ka(K, — Kj), 
Js =L,+ K3 Kz, 
Jo = 13+ Kglo+ K3(K,+ 53), 
Jp = Ly + K3Ly+ Ks (Ky + K,°), 
Ig =L,+ Kgly+ KgLy+ Ky KgLy+ (Ky+ K2) (Ky + Ky?) 


the invariant Ll, evidently being redundant. 


» 
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The converse process yields the formulae 


Tee 
oe eee ene 1 
Keele 
iene IR 


Ke=3l7.—1i+ 21) +16(, — 21,) +4 (R-341,4 2h), 
Wie ites tis halal 
L.= ita Ty , 
LJ, — Ia Js In (15 TDs 
L=J; Ih Js 1, (15 ts 
Jp PN Oy Fa) Iie 245) 
gy a ea) Poa) Sl Ta Dla la Da, 
Ee fp — 2h a 1a) (hh — 2 Je) 
ae Bl Ree dnd) A Lge tala) edges 
ee Cale Iola (lis = a) gg), 
In this case, to simplify the formulae resulting from the rearrangement of the 
symbols, it is necessary to make use of the Cayley Hamilton theorem (5.2), the 
formula obtained by putting e =a in (5.3), and the relation derived from this 
by interchange of a and b. Corresponding to the relations between J, and kK, 
in (7.7) and (7.8) there exist analogous relations between J} and K;. 
By eliminating L, from the expressions for J, and fj and substituting for 
K,, K,, L, (r=1 to 5, s=1 to 6) from (7.8), either J, or fg may be expressed 
as a polynomial in the remaining invariants (7.7). A straight-forward calculation 


ase A 2h= H+ @h— Mh ~~ 1) — 2h lo 
ik oh No Ae eRe 
eee ae 
ree err, ara] oe 
(21, + 1.215) + 

ere eRe ta 


( 
et) 
2 
et 
2 
io 
2 
1 
2 


8. A further reduction for transversely isotropic bodies 


By the procedure of the preceding section we may effect a considerable 
reduction in the number of coefficients and invariants of the type discussed in 
§ 5 which are necessary for the general representation of a transversely isotropic 
body. The steps in this reduction may be summarized as follows: 

(7) A given invariant J of the system of 3x3 matrices a, having the form 
(3.16) is expressed, if possible, as a polynomial in simpler invariants by means 
of reductions based on the Cayley Hamilton theorem and Peano’s theorem; an 
integrity basis for the isotropic invariants derived by this method is given by 
SPENCER & RIVLIN (1959a, b) and the additional invariants required for trans- 
versely isotropic bodies are given in § 5. 
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(ii) Each invariant of the 3 x3 system thus obtained is expressed by means 
of the symbolic method of § 7 as a polynomial in the invariants of the system 
(6.2) together with the leading terms of the constituent matrices, these leading 
terms being themselves invariant under the transversely isotropic group of 
transformations. 

(iii) The invariants of the system of 2x2 matrices (6.2) are expressed by 
means of the Cayley Hamilton theorem as polynomials in the invariants (6.6). 

(iv) By a reversal of step (7) again using symbolic methods, each invariant 
of the type (6.6) is replaced by a polynomial in the isotropic and transversely 
isotropic invariants of the system of 3 x3 matrices @,. 

(v) The resulting invariant system is simplified further, if possible, by applica- 
tions of the Cayley Hamilton theorem and its generalizations for a system of 
3 <3 matrices. 

This procedure yields a polynomial expression for J in terms of simpler 
invariants, provided that J is reducible. 

In each step of this reduction the equations employed are homogeneous 
throughout in the elements of the constituent matrices. In particular, in applying 
the process (77), if the given invariant of the 3 x3 system is of degrees /, m, n, ... 
in the matrices a, b, ¢, ... respectively, then each term of the relation connecting 
it with the invariants of the 2x2 system is also of degrees /, m,n, ... in the 
elements of a, b, c,.... A similar remark applies to the converse procedure (7v) ; 
the relations of § 7 illustrate this feature. Since each invariant of the 2x2 
system can be expressed as a polynomial in the invariants (6.6), and these in- 
variants are all of the fourth degree or less in the elements of the constituent 
matrices, it follows that any invariant of the 3 x3 system which is of degree 
greater than the fourth is reducible. Furthermore, each irreducible invariant of 
the 3 x3 system must be such that it arises as a result of applying the symbolic 
process to one or other of the invariants (6.6). Invariants involving two matrices 
have been examined in § 7. For the remaining invariants of (6.6) we have 


Ay Cap Og, = Ay 5 C55 gy — yy Cy 551 — Fy Cy + QO Cy, 
Aa xp Ap y by = 4; C;;45% Dyy— ayy C1 ;4;; by by; C4 j4— Ay 5 Cj14 4, 0%3+ (8.1) 
yg Oy Uy 5 Opt yy Ogg Cy g Apy tt Og yy My 5 Cp — M1 Oy Oy Ay - 
From these relations it follows that the only invariants of the 3 x3 system of 
matrices involving three or four different matrices which are not reducible are 
those of the types derived by taking the leading terms of (5.1) (I), (IJ), (III) 
and (5.17). If c,g=d,, the left hand member of the second of equations (8.1) 
_is reducible by means of the Cayley Hamilton theorem and this is sufficient to 


exclude the invariant formed by taking the leading term [a, a? a, |,, of the third 
product of (IIT) in (5.1). Thus 


[a be}, = yO; 5 Oj Cpr 
= Ay x Dap bp y lyr 1 Oy 4 Ox Cpr t C11 bp bpit 
+ Ay Om1 Oy p Spi + 41 O44 B14 la+ 44404 Oy gy Og + 


2 
+ 844643 Me Boa + 444 bri C44, 


(8.2) 
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and 
: Fy a Dag OB y Cy = B10 [Pay Opp + F5xy Op, bag — (dep)? ] Cy 


= 056 Oy Day Cyr + 31 061 [Opa Pag — (Bpp)?] - 
Again, invariants such as 


trabe, trabe?, trabed, 


which form part of the integrity basis for an isotropic system are all expressible 
in terms of simpler isotropic and transversely isotropic invariants by using 
stages (i), (2), (2v) of the reduction process described above. 

Summarizing the results of § 7 and the conclusions of the present section 
we find that any polynomial in the elements of the 7 symmetric 3 x3 matrices 
a@,,@,,...,@, which is invariant under the group of transformations (4.1) is 
expressible as a polynomial in the systems of invariants 


ia.) atta;. trae (i) 

tra;,a;, tr aja,, tra,a,  traja;, (ii) 

[@hx1, lags, (111) 

[a; Ghia, [a7 Gu, (a; GF) (2 +1), (iv) 

(a; 4; @,h11, (a; 4; 4, )11, [a; a; Qi li, (v) (8.4) 
[a; a, G;)11, [a, a, Gli, (1,7, & all different), 

[a;a; a, M11, [a; @;4,%,)1;, [4;4,4;Q]11, 

[a; a, @ Gli, (4,%4,4,)1, [4,4 Q)11, (vi) 


(1,7, k, 2 all different) , 


obtained by taking all possible choices of one, two, three or four of the letters 
7,7, k,l, as appropriate, from the numbers 1 to 7. In this scheme redundant 
products of type (v) are excluded by using (5.8). 

The invariants for transversely isotropic bodies are derived from (8.4) by 
replacing a, (p =1, 2, ..., w-+1) by the kinematic matrices E, A, (f =1, 2, ..., 1). 
The coefficients in the stress deformation relations follow by forming from (8.4) 
invariants which involve the transformation matrix C linearly. In this way we 
obtain the system of coefficients 


tO. te Cae tr as. (1) 

(Chi, [Ca@)ui, [Caihu, (11) 

(Ca; Gjji [Ca;ajju, [C a; Gj lh, [C a; Qi hy; (111) (8.5) 

[C a;a, Ahi, [Ca;,a,4;],, [Ca a; a, 45 

[Ca,a,a;],, [Ca,a;aji, [a,;Ca;a,)h1, (iv) 

(6 +7 +k+i) 

in which a,,a,;,q@, can represent any of the kinematic matrices E, A, (p=1, 
2,...,). The coefficients (i) in (8.5) are those appropriate to isotropic materials 
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and may be written in more familiar forms in terms of the convected system #. 
For we have 


oo” 08! ij 
ya eae 
1/00 a0 . 08 28'\—) _ (4 
te Oa — Free eat hoe ee Pa 


1/08 of an Of 
tO G5 = 5 laa “patches ral 
where Pye, (y'7 may represent the ae contravariant Pgs Suey 
from any of the kinematic tensors 7; ;, af’? with the metric tensor y’’, and 

represent the components of the corresponding tensor (that is, the elements iS: 
the matrix a,) referred to the cartesian reference frame x;. From (8.4) to (8. 6) 
it follows that the stress tensor t’’ for a transversely eae body may be 


written 
n+1 


49 Lf Ly i 1 1 1 2 
T= Ch Dot CHD [Pay OP +0, , a5, OP] + 
n+1 Fe ; ; ' @ 
ie) > [a are DY 1 Og, a, ay, BD}? a5 
Ges 1 m=1 
l~m 


n—1l n 


n+1 
ij 1 1 
ae Ci 2 ay 2 [' a a Ola. Dy mpt ( lane i i Dr pm aa 
=1 m=2 p=3 
l<m<p 


1 1 

1 ONG ae Ay Dytpt Pog ar a. ae Dipl ig (8.7) 
1 

at oe as, Oise Dy im] mG 


n+1 
os , i 
Sie Cy a > es Danas war Dimp si 


te cals as paseree fay 
GD cid, Uae SEES ate 


where Dy, D\”, Dl) D4, Dinp, %, Y%™ are polynomial invariant functions of 
the invariants derived from (8.4) by substituting the kinematic matrices E, A, 
for a,. By excluding from (8.7) the terms involving ®,, &, O*) , @,,,, and 
®,,,» and restricting %, Y% to be functions only of the invariants (i) and (11) 
in (8.4) we derive forms appropriate to isotropic bodies. These are not, however, 
the most general forms for the isotropic case, since by the procedure of §§ 7 and 8 
some of the isotropic coefficients and invariants have been expressed partly in 
terms of those peculiar to transversely isotropic materials. 


9. Elastic materials 


We consider an elastic material whose elastic properties are described in the 
usual way by means of a strain energy function W which is a polynomial in the 
strain components e,;. For such materials we have (GREEN & ZERNA (1954)) 

ot pea aw 
2VI 0e;; 06; 4 


I =|0;; + 2¢, ;|. (9.2) 


(9.1) 


where 
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From the formulation of §§ 4 and 5 we deduce that W can be expressed as a 
polynomial in the quantities 7, of (7.2) with e;; replacing a;;. Using (2.7) and 
(2.8) we then obtain 


Vin ae = ow 
v7 1) tf aj a nf! 
Trl BE 20 By t 3 ay 


(9.3) 
4 


- an ew ! ae an, ae a ow 
Ove On Ol, Coe Oe OM ONY Of, 


In the alternative formulation of § 6, the invariants K, of § 7 replace J, giving 
in place of (9.3) 


| 


T 


fit = ao aa ew (ee ay, at as ow , ae as) ow 


axe Ox OK, ' \Oxt dxP ° Oxb Oxe) ° OK, ' Ox Ox OK, 
an 3h | as api ow sy 
+ (Ge axe" xe sat) (eas onpe exe %1 5K) t (9.4) 
as 0) | as adi se a 
Ie  Ox% OxB oxB ca 1a 1B OK, é 


Equation (9.3) may evidently be regarded as a special case of (8.7), while (9.4) 
corresponds to the form which may be derived from (6.6). 


10. Orthotropic bodies 


In a previous paper (1958) the writer has derived for orthotropic bodies the 
relation 


1 = Aj} Oy + Brot Esp A She AR . pe : One OW. tf Lye AG (10.1) 


in which 
2 it 2 : 
Pe Ney Orgy yh sr Qi?) , (i) 
J- 1 4 sy 
Oe t= Ort ets oe}? , Qe. iS a 4 a x oe? , (11) (1 0.2) 
Ry 54 = Crs 19 OW + Cy: OW, (iii) 


54 18 equal to +14 or —1 according as 7, s,¢ is an even or odd permutation of 
4, 2,3 and is zero otherwise and in (10.1) summation is carried out over the 
repeated suffixes 7,s,¢. The functions 0;,, © are invariant functions of the 
quantities 


( 
Cris ay; , 
(1) ,(1) (1) 
Creer = UG Bi, =e, 
(1) (1) ,(1) (rh) (cL) (10.3) 
€12 &23 31> Cie R Aki» C7 Ak Ahi » A, 2 423 434 


(Gjek 1, 2)3) tay 4-4-1; %,7, & not summed), 


such that @{°) and @!§) change into each other by interchange of 7 and s 
(y--s-++t-r). The relation (10.1) is form invariant under the group of trans- 
formations 


Reet) = 1( x4, Xo; Bay (10.4) 
14* 
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describing the symmetry properties of orthotropic materials and (10.3) is the 
invariant system appropriate to the case where t'? depends upon the two kine- 
matic matrices e,,, a. 

The derivation of (10.1) to (10.3) depends upon the fact that in any product 
of the elements ¢;,, ayy which is form invariant under the transformation group 
(10.4), each of the suffixes 1, 2, 3 must occur an even number of times. If such 


products are arranged in the form 


b, LE eae Coe Co (1,7, k,l,...,m,” not summed), (10.5) 


aq Crk 


with equal suffixes adjacent, the limited range of values 1, 2, 3 for the indices 
i,],...,m, implies that any product (10.5) can be resolved into a product of 
the ivanont (10.3) in which at most three of the elements e; ;, a a‘) are multiplied 
together. An analogous argument limits the coefficients in (10.2). 

A similar procedure may be followed when 1’? depends upon any number of 
kinematic matrices ¢;,, ax’) (y =4,.2,...;”). In this case the invariant system 
corresponding to (10.3) may be derived from 


D1 9Co3 dey (10.6) 


by replacing each of b;;, c;;, d;; in all possible ways by e; ;, a av) (7 =1, 2,...,%), Ogee 
In place of (10.2) we obtain a set of functions 


Be ae 0, ? (1) 

One a by145 Di" ) (11) ) 
A (107, 

Re, == CnC =f Dare (On. = Crs) ? (iii) 


(P= 1523.03 Cees S862, HOw summed )= 


in which 0,,,c,, are replaced in all po ways by the kinematic matrices 
ey, 4) (p=1, 2,...,n) and OFF, OW" 7 are invariant functions of the 
invariants derived from (10.6) which are such that ¥4\"), y\) are interchanged by 
interchange of 7 and s. 

It should be noted that when any orthotropic invariant is expressed in terms 
of the basic system derived from (10.6), this resolution is not necessarily unique. 
For example, the invariant product 


I= By 9C3 431 frohoghsy, 
can be written 


I= (Or 2/2) (Cog ho3) - (dg, hy) 
= (01 2C93 453) - (fieMe3h31) = (fieCo34s1) - (01 2/93 Rg1) (10.8) 


re) 


and even if all of the constituent elements in J are equal, the first form of J in 
(10.8) represents a different resolution from the remainder. 
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Further Results in the Theory of Matrix Polynomials 


A. J.M. SPENCER & R. S. RIVLIN 


1. Introduction 
The theory of matrix polynomials and its application to the mechanics of 
isotropic continua has been discussed in earlier papers [J], [2], [3], [4]. The 
purpose of this paper is to give some further results which extend and in some 
respects simplify the theory which was put forward in [3] and [4]. 


The following theorem was proven in [3]: 


Theorem 1. Any matrix polynomial P in R 3X3 matrices Ap (P =1, 2, ..., R) 
can be expressed as a matrix polynomial of lower or equal extension and total degree, 
in which 

(1) each matrix product is either the unit matrix I or 1s formed from some or all 
of the factors dp (P =1, 2, ..., R) and at most two of the factors a} (P:=1,2; -- ee 


(it) no two factors in a single matrix product are the same; 


(111) each matrix product ts of lower or equal partial degree in each of the matrices 
Gp (P=1, 2,...,.R) than the matrix polynomial P; 


(iv) matrix products containing two of the factors a} (P =1, 2,..., R) contain 
them consecutively ; 


(v) no matrix product containing both of the factors ax and aj contains either 
of the factors Ax or az, unless it has the form ag az, az; 


(vt) Ax precedes az in any matrix product containing both ag and az. as factors. 


For definitions of the terms employed above and elsewhere in this paper, 
reference is made to [2] and |3]. 

The sets of matrix products which satisfy the conditions of Theorem 1 were 
enumerated in [3] for the cases R=1, 2,3, 4 and 5. Also in [3], a finite integrity 


basis for & symmetric 3 x3 matrices was determined, by a method based on 
Theorem 1. 


To this point all results were derived by means of algebraic operations based 
ultimately on the Hamilton-Cayley theorem. At this stage Peano’s theorem was 
introduced to further limit the number of elements in the integrity basis for R 
symmetric 3 x3 matrices, and also the number of necessary terms in the closed 
expression for a symmetric isotropic matrix polynomial in R symmetric 3 x3 
matrices, and the terms of this expression were enumerated. 
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In this note a stronger version of Theorem 1 will be proven, and this will 
be used to construct a closed form expression for a matrix polynomial P (not 
necessarily symmetric) in R 3 x3 matrices. The matrix products in this expression 
will all be of degree <5. For the case in which P is a symmetric matrix polynomial 
in R symmetric 3 x3 matrices, it will be shown that P can be expressed in a 
closed form which contains fewer terms than the corresponding closed expression 
which was given in [3]; moreover, all the results of this paper will be obtained 
only by algebraic operations based on the Hamilton-Cayley theorem, and without 
use of Peano’s theorem. 


Then, it will be shown that the integrity basis for R symmetric 3 x3 matrices, 
which was derived in [3] by the use of Peano’s theorem, and further discussed 
in [4] for the cases R<5, can be obtained directly from the results of this paper. 
Thus, it will be shown that all the results of [3] and [4], and some further results, 
are direct consequences of theorems of matrix algebra, and do not depend for 
their proof on general results of the theory of algebraic invariants. 


The application of these results to the formulation of constitutive equations 
describing the mechanical properties of isotropic materials for which the stress 
depends on a number of kinematic tensors is obvious. In § 11, they are applied 
to the case when the stress components are functionals of the displacement 
gradients, over the time interval from the instant at which the deformation 
commences to the instant of measurement of the stress. This problem has been 
discussed previously in [6]. 


2. Some results concerning matrix polynomials in 3 x 3 matrices 

In this section some results, proofs of which were given in [2] and [3], will 
be summarised for convenient reference. The notation Q(7) will be introduced 
to denote an arbitrary matrix polynomial of total degree less than or equal to 
y in the matrices ap (P =1, 2,..., R). If P is a matrix polynomial in these 
matrices, the statement P =Q/(r) is understood to mean “‘P can be expressed 
as a matrix polynomial of total degree less than or equal to 7 in the matrices @p 
(2402: >.., 2x) eeClearly: 

Q (7) Q(s) =Q(r +s). 


A similar notation was introduced by ADKINS [9]. 


In this notation, certain of the results given in [2] and [3] may be stated as 
follows. Let x, y and 2 be matrix products (not equal to the unit matrix I) of 
total degrees 7, s and ¢ respectively in the 3 x3 matrices ap (P =1, 2,..., R). Then 


xyztaucytyce+yer+eny+eyx=Q(r+s+i—1), (2.1) 
ay + yu?t+xyx=QO(27+s—1), (222) 

eye + xvyx = O(37+s—1), (2.3) 

xy 22? = O(27+s+ 2t—1), (2.4) 


Cy x2 = O37 25+ tA): (2.5) 
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In §§3, 4 and 5, we apply these formulas to show that matrix products of 
the forms 


(i) aga,ax,a,, gay aya,, ay AN ag; (2.6) 
(ii) pa, yGya%, Aga, ayazay, Aga, AS ay Ay, 

Ap A2a, Ay dy, Ag A, Ay Ay; (2.7) 
(iii) Ag A; Ay Ay As Ay (2.8) 


may be expressed as matrix polynomials of degree five or less in the matrices 
Ay, Az, Ay, Ay, As, Ap. 


3. Matrix products a, 4, @5,0\, Ac Qj AN,, Ai AX Ag Ay 
In (2.4), we replace a,y,2 by ax, @,, ax respectively and multiply the 
relation so obtained on the left by aj,;. This gives 


a, (Ag Ay Ay + Og AR a, + Gy Og Ay + ay, Ay Ag + ay ag a, +ay a, ax) =Q(5). (3.4) 


Now, by relations of the type (2.4), multiplied on the right by appropriate factors 
where necessary, each of the products a};a,aya@,, @),4, 4,0, @j,4xa,ay and 
ai,az,ax,a, can be expressed as a matrix polynomial of degree five or less in the 
matrices @x,@z,,@y and @y. Inserting these results into (3.1), it follows that 


yy Ay A a, + ah, ay ay ax = Q(5). (3.2) 


In a somewhat similar manner, replacing a, y and 2 in (2.1) by ag, aj,, and a} 
respectively, multiplying the relation so obtained on the right by a,, and using 
relations of the type (2.4) where possible, it is found that 


2 2 Pe 2 2 2 2 
Gy Ay Ax Ay + Ay Ay Aq Ay + Ag Ay Ay Ay + Ax Ay Ay a, = Q(5). (3.3) 
Next, in (2.1), replace # by a,ai,, y by ay, and 2 by ax, obtaining 


2 2 2 2 2 2 2 2 

Gy Ay Ay Ag + Ay Ay Ag Ay + Ay Ay Ay Ax + Ay Ag Ay Ay + (3.4) 

ee 2 2 : 
+ Ax Ay Ay Ay + Ax Ay A, Ay = Q(5). 

By relations of the type (2.4), multiplied on the left or right by appropriate 

factors where necessary, each of the matrix products a,aj,axay, @)Q,03,Ax, 

2 2 2 A 

ay 4x4, ay and axaya,ay,; can be expressed as a matrix polynomial of degree 

less than or equal to five in the matrices ax, @,,@,, and ay. Inserting these 

results in (3.4) gives the relation 


Ay Ay, AN Ag + Ag A, ay, AX = Q(5). (3-5) 


In a similar way, replacing #, y and @ in (2.1) by a}az, aj, and ag respectively, 
it may be shown that 


Ag Ay Ay A, + ay Ay a, ax = O(5). (3.6) 
From (3.3), (3.6), and the relation obtained from (3.6) by interchanging ay, and ay 


2 2: 2 2 2 2 
Ay Ay Ax Ay, — Ay Ay Ay, Ag + AY yy Ax Ay, — AX Ai, A, Ag = O (5). (3.7) 
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Hence, from (3.2), the relation obtained from (3.2) by interchanging ay, and ay, 
and (3.7), we obtain 


yy Ay Ax A, + AX ah, ay a, = Q(5). (3.8) 
Similarly, it may be-shown that 
Ax Ay Ay Ay + Ag A, AX a5, = O (5). (3.9) 
Next, in (2.5), replace x, y and 2 by ax+ay,, Gy and a, respectively. Then 
(Ax + Gy) Gy (Ax + ay)? a, = Q(5). 
Expanding, and using the relations (of type (2.5)) a,aj;a;a,=Q(5), and 
ay @y aya, =Q (5), this last relation becomes 
Og Ay Ag Ay Ay + Ay Aj, Ay Ay Ay, + Ay A5, Ag Ay A, + Ay a2, Ay Ag A+ 


(3.10) 
+} Og Qh AX Ay, + Ay Ay A a, = Q (5). 


Now, replace ay in (3.10) by —ay, and add the relation so obtained to (3.10). 
This gives 
Og Gy A, Ay + Ay C4, Ay Ay Ay + Ay G5, Ay Ag a, = O(5). (3.41) 


By relations readily derived from (2.2) and (2.4), and by (3.8) 


Ay Cay Ay Ay A, = — AX Ah, Ag Ap — Ay Ag AF, A, + O (5) =— a}, a5, ag a, +O (5), (3.12) 
and 
Ay Ais Ay Ag A, = — AY Ay Ag A, — Ay AN Ax A, + Q (5) =Q (5). (3.13) 


Inserting (3.12) and (3.13) in (3.14), the relation 

Og Dy Ay Ay, — AY Ay Ax A, = Q (5) (3.14) 
is obtained. Hence, from (3.6) and (3.14) 

yy AX Ay A + AY A}, Ag A, = Q (5). (3.15) 
However, from (3.2) and (3.8), 

a2, a>, A, Ay — Ay ai, Ag A, = (5), (3.16) 


and it therefore follows that 
hy Ay A, Ax = Q(5). (3.17) 


From (3.17), (3.5), (3-6), and relations obtained from these by permuting the 
matrices, it is readily shown that 
ay aya d,=O(5), ax ayaya,=Q(5),  axapayay =Q(5). (3.18) 
We thus have 


2 22 
Lemma 1. Any matrix product of the form @yaxaga,, Anaya, or 
A A, Ai, A 1s expressible as a matrix polynomial of total degree five im the matrices 


Ag, A,, Ay, Ay. 
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4, Matrix products a, @y Gy @%, Ag, Ay A3y, Ap A, AZ Ay, Ay, 
Ac AEA, Ay Ay, U3 Ag Ay Ay Ay 
Interchanging a,, and @y in the first of the relations (3.18) and adding the 
relation so obtained to (3.18), we obtain 
(ah, @y + Gy Gy) Ax a, = Q(5)- 
Using a relation of the type (2.2), this leads to 
Ay Ay Ay Ax A, = Q (5). 


Applying similar considerations to the second and third of the relations (3.18), 
the three relations (3.18) yield 


Ay, Ay Ay Ax a, = Q(5), Ag Ay, Ay Ay A, = Q(5), Ay A, Ay Ay Ay = Q(5). (4.1) 


Also, from relations of the type (2.4), 
Ai Ax Aya, =Q(5), Ai axa,ay=Q(5), ax aya,ay =Q(5). (4.2) 
In the first of the relations (3.18) replace ay, by as and ay by @y,+ay. Then 
ds (yz + Ay Ay + Ay Ay + AN) Axa, = Q(5). (4.3) 


By relations of the type (3.18), ataj;a,a@,—=@Q(5) and atayaza,=Q (5). Hence 
(4.3) yields 
aS (Ay Ay + Ay Ay) Ax a, = Q(5). (4.4) 


Alternatively, replacing @y by @,,+a,, and ay by dg in the first of (3.18), and 
proceeding in a manner similar to that used in deriving (4.4), the relation 
(yy + Ay Ay) AS axa, = Q(5), (4.5) 
may be obtained. In a similar manner, we obtain from the second and third 
of the relations (3.18) the four additional relations 
Ag US (Ay Ay + Ay Ay) @, =Q(5), Ax (Qy, Ay + Ay Ay) aa, = Q(5), 


Ay A, aS (yy, @y + Ay Gy) =Q(5), Ax A, (Ay, Ay + Ay Ay) az = Q (5). 


(4.6) 


Again, replacing @y by @y+ dy and ay by ag in (4.1) and employing relations, 

of the types (2.2) and (3.18), we obtain 
(Ay AS dy + Ay AS Ay) Ax @, = Q(5), 

Rt Ax (Ay AS Ay + Ay AS Ay) a, = O (5) (4.7) 
Ax A, (Ay AS dy + Ay AS ay) = Q(5). 
Replacing ay by as and ay by ay,+ay in (4.2) and employing relations of the 
types (4.2), we obtain 
AS Ax (Ay, Ay + Ay Ay) a, = Q(5), 
AS Ag Ay (Ay Ay + Ay Ay) = O (5) (4.8) 


and 2 
A, Ws My (Ay, Ay + Ay Ay) = O(5). 
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Similarly replacing ay by a; and ay by ay+4ay in (4.2), we obtain, with (4.2), 


(@y Ay + dy ay) ax asa, = Q(5), 
(@y Ay + Ay Ay) Ax a, aS = Q (5) (4.9) 


and Ax (Ay Ay + Ay Ay) aa = Q(5). 


The fifteen relations given in equations (4.4) to (4.9) lead immediately to 
the result that each of the five matrix products a,a@,ayaya%, ..., a2ag4, Ay ay 
is equal to the sum of the negative of the product formed from it by interchanging 
a pair of adjacent matrices in the sequence ax, @,, @y, @y and a matrix poly- 
nomial of total degree five or less. From this result we obtain immediately 


Lemma 2. Any matrix product of 3X3 matrices dp (P =1,2,..., R) which 
has the form xaty, where xy is a matrix product of the form Aga, AyaAy, is equal 
to any matrix product formed from xazy by an even permutation of matrices in 
the sequence Ax, A; , Ay, Ay and to the negative of any matrix product formed from 
xazy by an odd permutation of matrices in the sequence Ag, A, Ay, Ay, together, 
in each case, with a matrix polynomial of total degree five or less in the matrices Ap. 

It will now be shown that each of the matrix products aga, ayQya%, ..., 
a$a,4,@y4y can be expressed as a matrix polynomial of total degree five or 
less in the matrices ap. In (2.4), replace x, y, 2 by a;, €x, Ay+ Ay respectively. 
We obtain 

AS A (iy + Ay Ay + Ay Ay + ay) = Q (4). 


Since, by relations of the type (2.4), aga,aj,=Q(4) and afa,ay =Q (4), this 
may be written 
AS Ax (Ay Ay + Ay Ay) = Q (A). (4.10) 


In (4.10), replace @x by aga, 

AS Ag My (Ay Ay + Ay Ay) = Q (5). (4.14) 
Also in (4.10), replace ay by a, ay, 

AS Ax (A, Ay Ay + Ay a, Ay) = Q (5). (4.12) 


Again in (4.10), replace ay, by a,, and multiply the relation so obtained on the 
right by ayy 


AS Ax (Az Ay + Ay G;) Ay = Q(5). (4.13) 
Adding (4.11) to (4.12), and subtracting (4.13), the relation 
AS Ag A, Ay Ay = Q (5) (4.14) 
is obtained. In a similar way, it may be shown that 
A Ay Ay Ay as = Q(5). (4.15) 


Now, in (2.1), replace # by a&, y by axa, and 2 by a),ay, obtaining 
2 2 
a2 Ag Gy, Ay Ay + Az Ay Ay Ag Gy, + Ag A, AS Ay, Ay + Ag A, Ay Ay $+ 


+ dy Ay A Ag A, + Ay Ay Ax GAS = Q(5). 


220 A. J. M. SPENCER & R. S. RIVLIN: 


Inserting relations of the types (4.14) and (4.15), this becomes 
Oz, A Ay Ay + Ay Ay 5a, a, = O(5). (4.16) 
Since @,,@ya2a,a,, can be obtained from a, @,a$ay@y by an even permutation 


of the matrices a, @,, @, and dy in the product @,¢4@,4),;@y, it follows from 


(4.16) and Lemma 2 that Og Oy, Fy, yy — O15), (4.17) 


Next, in (2.1), replace a by a%, y by @ga,@y, and z by ay. Of the six pro- 
ducts which appear in the resulting relation, four can be expressed as matrix 
polynomials of lower degree by means of relations of the types (4.14) and (4.15). 
The relation then becomes 


Ay Uy Ay AS Ay + Ay AS Ag A, Ay, = Q (5). 
Applying Lemma 2, this relation may be written 
A Ay My AS Ay — Ax AS Ay Ay, Ay = Q (5). (4.18) 
Again in (2.1), replace « by a%, y by axa@,, 2 by @),, and multiply the resulting 
relation on the right by ay. We obtain 
(a3 ag Ay Ay + AS Ay Ag A, + Ag M, AS Ay + Ag M, My AS + 
+ @y 3 AG, + Ay Ax d, AS) dy = Q(5). 


Inserting relations of the types (4.14) and (4.17), we have 


2 2 2 
Ag A, Ay AS Ay + Ay, AS Ag A, Ay + Ay Ag , AS Ay = Q (5). 


With Lemma 2, this last relation may be written 


20g Ay Ay, US Ay + Ax ASM, Ay Ay = Q(5). (4.19) 
Hence, from (4.18) and (4.19) 
Ac ASA, Ay Ay =Q(5), Ag A, Ay, AS ay = Q(5). (4.20) 
Equations (4.14), (4.15), (4.17) and (4.20) express 
Lemma 3. Any matrix product of 3X3 matrices Ap (P =1, 2,..., R) which 


2 
has the form xasSY, where x“y isa matrix product of the form Axa, ayy, is equal 
to a matrix polynonial of total degree five or less. 


5. Matrix product a, @, @y, Ay a, ay 
In (4.14), replace a, by ag+-a;. We obtain 
(a5 + Gs yp + Ap As + AF) Ag A, Ay Ay = Q(5). 

Employing relations of the type (4.14), it follows that 

(5 Ap + Ap As) Ax G Ay Ay = Q(5). (5.1) 
In a similar manner, each of the relations (4.15), (4.17) and (4.20) leads to a 
relation in which a% in the original relation is replaced by a,a;-+-ara,. Thus, 
AM, Ay Ay Asay is equal to the sum of the negative of the matrix product formed 


from it by interchange of any two adjacent factors and a matrix polynomial of 
degree five or less. We obtain immediately 
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Lemma 4. The matrix product Ag, Ay Aya, formed from six 3 X3 matrices 
is equal to any matrix product formed from it by an even permutation of its factors 
and to the negative of any matrix product formed from it by an odd permutation 
of wuts factors, together, in each case, with a matrix polynomial of degree five or less 
in the matrices. 

To complete the discussion of this matrix product, in (2.1) replace x, y and z 
by @g@,, @y,ay and aay respectively. This gives 


Ox A, Ay Ay As Ap + Ag Ay As Az Ay Ay + Ay Ay Ag A, As Ap + Ay Ay As Ap Ax A, + 


+ As Ap Ax A, Ay Ay + As Ap Ay Ay Ax A, =Q (5). 


(5.2) 


It is readily verified that each of the matrix products on the left hand side of 
(5.2) can be obtained from a@,a,ay@ya,a; by an even permutation of the 
factors. Hence, from Lemma 4 and (5.2) 


Ag M, Ay Ay Asa, = Q(5). (5.3) 
This expresses 
Lemma 5. The matrix product AQ, Ay Ay asap formed from six 3 x3 matrices 
1s equal to a matrix polynomial of degree five or less. 


We note that Lemma 3 and Lemma1 may be obtained from Lemma 5 by 
taking appropriate pairs of the factors in the matrix product to be equal. 


From Theorem 1 and Lemmas 1, 3 and 5, we obtain immediately 

Lemma 6. Any matrix product of degree six or more in the 3X3 matrices Ap 
(P=1,2,..., R) may be expressed as a matrix polynomial of degree five or less 
in these matrices. 

An alternative derivation of Lemma 6 is given in §12. Although this is 
shorter than the derivation given above, it fits into the general framework of 
the discussion somewhat less well. 


6. Some matrix products of degree five 


In this section we shall consider those matrix products which satisfy the 
conditions of Theorem 1, contain a factor ag and a factor a% and have total 
degree five or less. These have the forms 


; 2 ay, , Ag Ay, A (6.4) 
Qn, 4K, Ax, Akay, Axa, ay Ax, Lan Ay aK. 
In a manner similar to that in which (4.10) was derived, we obtain * 
2 
(@y, Ay + Ay Ay) Ag as = Q(4). (6.2) 


In (6.2), replace ay, @y, dx and a, by Gx, ,, @y and x respectively, obtaining 


the relation 
2 2 A 
Ag A, Ay Ak + A, Ag Ay Ax = Q(4). (6.3) 


Also, by a relation of the type (2.3) 


Ax A, Ak Ay = — AK A, Ax Ay + Q(4), (6.4) 


* Relations equivalent to (6.3) and (6.4) have been given by ADKINs [4] 
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and by a relation of the type (4.10) 

Oi A, Ax Ay = — Ag My, Ay Ax + QA). (6.5) 
Further, by a relation similar to (2.3) 

Oj A, Ay Ag = — Ag A, Ay Ae + Q(4), (6.6) 
and thus, from (6.4), (6.5) and (6.6) 

Ay 0, 0, Ay = — Ag Uy, Ay a + O(4). (6.7) 


Equations (6.3) and (6.7) show that each of the matrix Abia? (6.1) a 
be expressed as the sum of one or other of the products AA, Ax and Ag A, Ay, Ax 
(with possibly the coefficient —1) and a matrix polynomial of lower total degree 
than that of the original matrix product. 


Combining this result with Theorem 1 and Lemma 6, we obtain 


Theorem 1’. Any matrix polynomial in R, 3 x3 matrices Ap (P =1, 2, ..., R) 
may be expressed as a matrix polynomial of lower or equal extension and total degree, 
in which each matrix product satisfies the conditions (1) to (v) of Theorem 1 and in 
addition 

(vi) each matrix product which contains both ax and ak as factors has Ag as 
the first factor and ai as the last factor; 


(vit) no matrix product has total degree greater than five in the matrices Ap 
(W278 ee eae 39 
7. Matrix polynomials in R 3x 3 matrices 
In this section Theorem 1’ will be applied to obtain an explicit closed form 
expression for an arbitrary matrix polynomial in & 3 x3 matrices. Enumerating 
all the matrix products which satisfy the conditions (¢)—(v), (vz’) and (viz) of 
Theorems 1 and 1’, it is seen that any matrix polynomial P in R, 3 x3 matrices 
ap (P =1, 2,...,.R) can be expressed as a matrix polynomial 


P=) 41 
where ¥; are scalar polynomial invariants of the matrices @p, and II; are matrix 
products of the forms 


; 2. 2 2 29 
T; Gg, Qk; Axa,, axay, aka, axaz, Ay My, A, Ax A, Ax; 


2 2 gen” Abo “s 2 2 a; 
Ux Ady, Aka, Ay, Aq Apay, And, Ay, AcAL Ay, Axa, Ay, Ax A, Ay Ax; 


(7.1) 


2 3 2 
An A, Ay Ay, AKA, Aydy, AA; AyAy, Axa, AiAy, Aga, ay az; 


Ax Ay Ay Ay As ; 


and where K, L, M, N and S are integers, which are all different, and are chosen 
in all possible ways from the integers 1, 2,..., R. 

The sets of matrix products in (7.1) are included in the sets of matrix products 
which satisfy the conditions of Theorem 1, for R<5, and which were enumerated 
in [3]. For R=1 and for k =2, the set of matrix products given by (7.1) is the 
same as the corresponding set given in [3]. For R=3, 4 or 5, the set given 
by (7.1) contains fewer members than the corresponding set given in [3]. 
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It may be noted that certain of the matrix products of the forms (7.4) can ~ 

be expressed in terms of other matrix products of these forms. The arbitrary 

matrix polynomial P in R, 3 x3 matrices ap may therefore be expressed as a 

matrix polynomial involving fewer matrix products than are implied by (7.1). 

We shall, however, leave the discussion of such further limitations on the form 
of the matrix polynomial to a later paper. 


8. Symmetric matrix polynomials in R symmetric 3 x 3 matrices 


Suppose now that P is a symmetric matrix polynomial in symmetric 3 x3 
matrices: @p (P =41,2,.,.,R). Then 


P=3(P+P), 
where P’ is the transpose of P. Hence, by the results of § 7, P can be expressed 
in the form ; 
P=3)>4,(0,+1), (8.1) 


where y; and II; are defined as in § 7, and II; is the transpose of II;. Since the 
matrices @p are symmetric matrices, II; is derived from II; by reversing the 
order of the factors in II,. Now, it follows from relations of the type (2.3) that 
2 2 
Qn a, Ak + Ae a, Ag = . 
and Sa aS a (8.2) 
Ak A, Ak + Ax a, Ax = Q(4). 


From § 7, together with (8.1) and (8.2) it follows that P may be expressed as a 
sum, with scalar polynomial invariant coefficients, of terms of the forms 


T; az, ay; aga,+apaz, ataz;ta,az, aaj -+azaz; 

Ax A, Ay + AyA,dgn, AA, Ay + Ay, a, Ag A, Ay + Ay, Aj Ax, 
AKA, Ay + Ay A, AK, AA, Ay Ae + Akay A, ax; — (8.3) 

Ax A, Ay Ay+AyAy A, Ag, AKA, Ay Ay+Ay Ay A, AK, Ax A, Ay Ay + Ay Ay, A; Ax; 

Ax A, Ay Ay As + As Ay Ay A, Ax, 


for all selections of K,L, M,N and S (all different) chosen from the integers 
12 ae ee Ae 

In [3] it was proven that a symmetric matrix polynomial P in symmetric 
3X3 matrices @p may be expressed as a sum, with scalar coefficients, of terms 
which were enumerated in Theorem 5 of [3]. The set of terms given by (8.3) 
is included in the set which was given in [3], but contains fewer members; the 
omission of certain terms which appear in Theorem 5 of [3] but not in (8.3) 
above is a consequence of (vz’) and (viz) of Theorem 1’ of this paper. Moreover, 
the set (8.3) has been derived from the Hamilton-Cayley theorem by methods 
of matrix algebra, and does not require the use of Peano’s theorem, upon which 
Theorem 5 of [3] is partly based, for its derivation. 


9. Invariants of R symmetric 3 x 3 matrices 
It is easily verified that all of the results regarding the invariants under the 
orthogonal transformation group of R symmetric 3 x3 matrices which were 
obtained in [3] and [4] can be obtained from Theorem 1’ by a slight modification 
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of the arguments used in [3] and [4]. As shown in [3], any polynomial scalar 
invariant of the R symmetric 3 x3 matrices ap (P =1, 2, ..., R) may be expressed 
as a polynomial in invariants of the form tray I;, where II, +I; are the expres- 
sions (8.3). II; is of total degree less than or equal to five in the ap and it follows 
that tra,, II, is of total degree less than or equal to six in the ap. Thus every 
irreducible invariant of the symmetric 3 x3 matrices @p is of degree less than or 
equal to six in the ap. 

Also, any invariant of degree six in the ap either involves five or fewer of 
the matrices, is reducible, or is equivalent * to some linear combination of the 
invariants tr @¢@,@),@ya,a,. These results are sharper than those obtained in 
[3], where it was shown by an argument based on Peano’s theorem that every 
irreducible invariant of the ap is of degree less than or equal to seven in the ap 
and any invariant of degree seven either involves six or fewer of the matrices, 
is reducible or is equivalent to some linear combination of the invariants 
tr dy A; Ay Aya, Ap ay. 

The problem of finding an integrity basis for any number of symmetric 3 x3 
matrices is thus reduced to that of finding an integrity basis for five symmetric 
3 x3 matrices. 


10. Integrity basis for five symmetric 3 x 3 matrices 

In a previous paper [3], finite integrity bases for five or fewer symmetric 
3x3 matrices, under the orthogonal group of transformations, were derived. 
In [4], it was shown that certain elements of these bases are reducible, and 
equivalences were established between other elements of the bases. Integrity 
bases containing fewer elements were thus obtained. In particular, it was shown 
that an integrity basis for five symmetric 3 x3 matrices a, b, c¢, d and e consists 
of the integrity bases for the five sets of four matrices which can be selected 
from a, b, c,d, e, together with the invariants 


(i) trabede, trabdec,  trabeecd, 
tracd be, tracbed, trad bee; 


(ii) trabcede*®, treabde’?, tredabe?, 
tracdbe’®, trdacbe*®, tradbce® 


and invariants obtained from these by permuting a, b, ¢, d, e cyclically ; 


(iii) trabecde’?, trbceade®, treaebde?, 
trbcedeae*®, trebaede®, treadebe?, tredbeae 


and invariants obtained from these by permuting 4, b, ¢, d, e cyclically; 


(iv) trabcd@e®?, trabde’e, tradeb?e?, trdbeca e?, 
trabee’d, traceb’d’, trebed a?, 
tradeb’c?, trbde@e, tredea?b?. 


In this section it will be shown that the number of elements in the integrity 
basis can be further reduced. 


; e An invariant J is said to be reducible, if it can be expressed as a polynomial 
in invariants of degree less than the degree of I and this property of J is denoted 
I=0. If J, and J, are invariants and J,— J,=0, then I, and I, are said to be equi- 
valent, and this is denoted J,=/,. 
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The following results, proofs of which were given in [3], will be used. Let 
a, b,c, x and y by 3 x3 matrices, not equal to the unit matrix J. Then 


try(abc+acb-+bac+bea+cab+ecba)=o0, (10.1) 

png, try@’azb?=0. (10.2) 
Also, if a, b and € are symmetric 3 x3 matrices 

tra’b?e?=0. (10.3) 


Use will also be made of the following Lemmas, which are consequences of the 
definition of matrix multiplication: 


Lemma 7. The trace of a matrix product formed from 3 X3 matrices 1s unaltered 
by cyclic permutation of the factors in the product. 


Lemma 8. The trace of a matrix product formed from symmetric 3 <3 matrices 
1s unaltered if the order of the factors in the product is reversed. 

It follows from § 9 that the trace of any product formed from the five sym- 
metric 3 x3 matrices a, b, c, d, e is reducible, if it has degree greater than six. 
Since each of the invariants (iii) and (iv) has degree seven, they are all reducible 
and may be omitted from the integrity basis. 


We now consider the invariant 
trabede’ (10.4) 


and invariants formed from it by permutation of a, b, c,d and e. The set (ii) 
consists of selected invariants of this type. It was shown in [4] that any invariant 
of this type can be expressed in terms of those given in (11). We shall now derive 
a number of equivalence relations between these invariants, which enables us to 
construct an integrity basis involving fewer invariants of the type (10.4) than 
are included in (ii). 

First, suppose that a, b, c, d and e are symmetric 3 x3 matrices and in (10.2) 
replace y, a, x and b by a, b+ e, d and e€ respectively, obtaining 

tra(b?+be+cb+ec’\de=0. 
Since, by relations of the type (10.2), tra b’de?=0, and tr ac’de?=0, this last 
1 b itt 

relation may be written Gath mebyd @=0: (10.5) 


In (10.3), replace a, b and ec by a+b, c+d and e respectively. Then 


tr(a+ b)?(e+ d)?e=0. 


Expanding the above relation, and inserting relations of the type (10.3) wherever 
possible, this gives t+ (ab + ba) (ed + de) e=0. (10.6) 

Relations between invariants derived from (10.4) by permutation of a, b, c,d 
and e may also be obtained by making various substitutions in the relation (10.1). 
It will be shown that any relation obtained in this way can also be derived from 
relations of the types (10.5) and (10.6). For, setting y=de? in (10.1), and 
grouping the factors in the indicated manner, we obtain 


trdea(be+cb)+trde’b(ac+ca)+trdec(ab+ba)=0, (10.7) 
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a relation which clearly can be obtained by forming the sum of three relations 
of the type (10.5), and applying Lemma 7. Further, replacing y, a, b and ¢€ in 
(10.1) by ab, c,d and e? respectively, the relation 


trab(edet+ced+dce’+dec}+ecd+edc) === ())5 (10.8) 
is obtained. Using Lemmas 7 and 8, (10.8) may be written 
tr(ab+ba)(cd+dec)e?+tre(ab+bajde’=0, (10.9) 


and this result may also be obtained by adding (10.6) to a relation of the type 
(10.5). It may be verified that any relation between invariants of the type (10.4) 
which can be obtained by making substitutions in (10.1) may, by use of Lemmas 7 
and 8, be expressed in one or other of the forms (10.7) or (10.8), or a form obtained 
from (10.7) or (10.8) by permuting the matrices a, b,e,d and e. Hence the 
relation (10.1) can yield no equivalences between invariants of the type (10.4) 
which are not contained in relations of the types (10.5) and (10.6). In considering 
equivalences between invariants of the form (10.4), therefore, only (10.5) and 
(10.6) and relations obtained from them by permuting the matrices will be used. 

There are six distinct relations of the type (10.5) which are obtained by 
permuting a, b, ce and d in (10.5). They are 


trabcede?+tracbde=0, trabdce?+tradbce’?=0, 
tracdbe?+tradcbe’=0, treabde?+itrcbade?=0, (10.10) 
trbacde?+trbcade=0, trbadce?+trbdace’=0. 


Also, there are three distinct relations of the type (10.6) which are obtained by 
permuting a, b, c and d in (10.6). These are 


trabcde?+trabdce?+trbacde?+trbadce’=0O, (10.11) 
tracbde’?+tracdbe?4+ treabde?+trbdace?=0, (10.12) 
tradbce®?+tradcbe?+trebade®?+trbcade?=0. (10.13) 


In (10.12) and (10.13) the order of the factors in some of the invariants has been 
altered in accordance with Lemmas 7 and 8. 

It is seen that the sum of the terms on the left of the six relations (10.10) 
is equal to the sum of the terms on the left of the three relations (10.11), (10.12), 
(10.13). Hence, not all of the nine equations (10.10) to (10.13) are independent, 
and one of them, say (10.13), can be regarded as a consequence of the remaining 
eight equations, and need not be considered further. By considering the matrix 
of the coefficients of the twelve invariants occurring in the eight relations (10.10), 
(10.11) and (10.12), it is found that these equations are independent, and may 
be used to express eight of the invariants, suitably chosen, in terms of the other 
four invariants. For example, (10.10), (10.41) and (10.12) may be used to express 
tracbde’, tradbce?, tr acdbe’, tr cabde?, tr ebade?, tr beade?, tr badce 
and tr bdace? in terms of trabede?, trabdce?, tradcbe? and tr bacde?. 
Hence of the invariant (10.4) and invariants derived from (10.4) by permuting 
a, b,c and d, it is necessary to retain in the integrity basis only these last four 
invariants. Therefore, of the invariant (10.4) and invariants derived from (10.4) 
by permutation of the matrices a, b, c,d and e, it is sufficient to retain in the 
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integrity basis only the invariants 
(ii’) trabcde’?, trabdce?, tradcbe?, trbacde? 


and invariants obtained from these by permuting the matrices a, b, c,d and e 
cyclically. We thus have the conclusion that an integrity basis for the five sym- 
metric 3 X3 matrices a, b, c,d and e may be formed by the integrity bases for 
the five sets of four matrices which can be selected from a, b, ce, d and e, together 
with the invariants enumerated in (i) and (ii’) and, invariants obtained from 
(ii’) by permuting a, b, c,d and e cyclically. 


11. Application to the mechanics of materials with memory 


In a previous paper [6], GREEN, RIVLIN & SPENCER have discussed the form 
taken by the constitutive equation for deformations of an isotropic material if 
the stress components in an element of the material at the instant of time f, 
referred to a rectangular Cartesian coordinate system %, are assumed to have 
continuous functional dependence on the deformation gradients in it at all times 
up to the time ¢ and particular dependence on the deformation gradients at the 
instant ¢. The deformation is described by the dependence of the coordinates 
x,(t) of a generic particle of the body on the time 7 and its coordinates X,; when 
t=0. 

The following notation is employed: 
ae ele ee Le 
| g=g), ¢=[6x/6X;|), (11.1) 


c’ = transpose of ¢. 


4, = % 


g(t) = lle; (7) 


With this notation it was shown that the stress matrix J’ in the coordinate 
system % is given by 


T2686) (14.2) 

where 
6 tt t 
F epg WS se, J i) Per Ona Coes: Ti) (Xie 4) d%4%...dT,. (11.3) 
v v=1 » 0 0 0 

In (11.3), x” for the various values of y are the matrix products in Table 1 which 
have all of the matrices g(t,), g(t), -.-,g(t,) as factors. x!” is the transpose 
of x. 3,(t, T:, -.., T,) is a function of ¢, %, T2,..., T, and a polynomial in the 
invariants 


tt t 
tr @®) and eles eerste) Os eee tae (Bois 57) © (44) 


in which tr @ are the invariants in Table 2 which involve none of the matrices 
g(t) (8 =1, 2,..., 7) and tr @{” are the invariants in Table 2 which involve all 
of the matrices g(t), J (T2),-.-, J (T,)- 

By an argument similar to that employed in deriving this result and using 
the results developed in previous sections of the present paper, it follows that 


the stress matrix J’ may be expressed in the form (11.2), where F is given by 


5 tat t 
F= D9, ~O 420 0 Sf L(t, te, es Ty) P+ Katya ty...dt,. (11.5) 
v p=1 7 00 0 


15 
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x” are now the matrix products above the line in Table 4 which have all of the 
matrices g(t), 9(t2),---, 9(t,) as factors. 8,(t%, T2, --+5 Ty) is a function of ¢, 
Utes rey oe a polynomial in the invariants (11.4) with 8 =1, 2,...,6. We 
note that ee 6°) (B =0, 1, ..., 6) are now the invariants above the line in east: ps 


Table 1 


( 

99 (%) G(T); 9 (%) 99 (2), 979 (%) G(T), 9(%) 9?G (tT), 
( 
( 


9(%), 9(%)99(T2) G(T), G?G(%) G(T) G(T), 


9 (%) 979 (72) J (Ts) ; 
99(%) G(T) 9(%) G(T), G(H)IG (Tz) 9 (T3) G(T), 
9 (%%) 9 (T2) Gg (Ts) J (Ts) 
9(%), 9(%)G(%), G9(%) G(T) G(T), G(t) G(T) 9 (Ts) Y (Ta), 
9 (%) 9 (2) 9 (Ts) G (Ta) Y (Ts) 
9 9(%) 97g (T2); 
99 (%) 9 (7) 9 (Ts) 9”, 9 (%) 99 (T) g (Ts) g*, 
9 (7%) 9 (72) 9G (Ts) 9”, 9 (7%) 99 (T2) 9°G (Ts); 
99 (7%) G(T) G(T) G(T), G(T) 9?G(T2) G(T) F(T), 
9 (%) G(T) 97g (Ts) G (Ta); 
99 (7%) G(T) (Ts) G(T) G (Ts), G(T) PG (T2) Y (Ts) Y (Ts) G (Ts) 
9 (%) 9 (T2) FG (Ts) G(T) (Ts) 5 
9 (7%) 9 (T2) F(T) G (Ts) Y (Ts) J (To) 
Table 2 


trg(t), trg(%)g(t), trg(t%) g(t.) G(T), 
T3)9(T4), tr g(t) G(T2) G(T) G (Ts) Y (Ts), 
) 9 (Ts) F(T) G(T) 5 


a 


1) 9% trg (ty) 9(t2) 99(%) g(t) 9°, 


12. Appendix. An alternative derivation of Lemma 6 
Lemma 6, which expresses the fact that any matrix product of degree six 
in the @p may be expressed as a matrix polynomial of degree five or less in the 
ap, may also be derived in the following way. It was shown in [3] that any 
invariant of the form tr ata, a%,y is reducible; that is, in the notation of [3] 
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and [4 

[4] tra,4,a7 y= 0, (y +1). (12.1) 
Replacing @y by @y-+ ay, and using the relations tr aja, a3;y =0, tr aza,ary 
=0, this becomes > 
tr a§ x (Ay Ay + Gy ay) y = 0. (12.2) 
Now in (12.2) replace @x by aga,. This gives 


tr @§ Ax @; (Ay, Ay + Ay Ay) y = 0. (12.3) 
Next, in (12.2), replace ay by a, ay, 
tr a2 dg (a, Ay dy + Aya, ay) y = 0. (12.4) 
Also in (12.2), replace ay by a,, and y by ayy 
tr a3 dg (a, dy + Ay a,) ayy = 0. (12.5) 
Adding (12.3) to (42.4), and subtracting (12.5), it follows that 
tr az dy a, ay ayy =0. (12.6) 


Since the trace of a matrix product is unaltered by a cyclic permutation of the 
factors of the product, it follows by making various substitutions in (12.6) that 
trad, ata,ayayy=0, trad, atayayy =0, 42.7) 
2. 
trad, Ay asayy=0, traga,ayayaty =0. 
In for example (12.6), replace a; by a,;+a,y, and insert the relations 
tr aad, ay ayy =0, tr apaga,ayayy=0. This gives 


tr (ds Ap + Az a;) Ag A, Ay ayy =0. (12.8) 


In a similar way, each of the relations (12.7) leads to a relation in which a% in 
the original relation is replaced by a,a;+arya,. It follows that any invariant of 
the form trax @,@y Ay 4,@7y is equivalent to the negative of any invariant formed 
from tr d@g,@,@yaya,a;y by the interchange of any two adjacent factors from 
Ay, Az, 4, 4y, 4, and ar. Hence any invariant of the form tr a, @,@y@ya,ary 
is equivalent to any invariant formed from it by an even permutation of the 
factors Gx, A, Gy, Ay, and a;, and to the negative of any invariant formed 
from it by an odd permutation of these factors. 

It was also shown in [3] that if a, b,c and y are 3 X3 matrices, not equal 

I, th 
Ont auret tr(abe+acb+bac+bea+cab+cbajy=0. (12.9) 
In (12.9), replace a, b and € by aga,, @y@y and a,ar respectively, obtaining 
tr (Ax Ay Ay Ay As Ap + Ax Ay As Ap Ay Ay + Ay Ay A A, As Apr (12.10) 
+ Ay Ay As Ap Ag A, + As Ap Ay Az Ay, Ay + As Ap Ay Ay Ax d,) Y = 0. 


It may be verified that each of the invariants in (12.10) may be derived from 
tr Ag Q,AyAyadsary by an even permutation of @,,4,,@y,@y,d, and ay. 
Hence each of these invariants is equivalent to tr dg @,Qy@y4,a,y, and it follows 


from (12.10) that tr Og Ay Ay Ay Ag py = 0. (12.11) 
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Now denote the invariant tra,@,@ya@ya,ary by J. By (12.11), J is a 
reducible invariant, and so can be expressed as a polynomial in invariants of 
lower degree. From the manner in which (12.11) is derived, it follows that the 
expression for J must be of the form 


J =X GaSe: (12.12) 


where the g, are invariants which are polynomials in traces of matrix products 
of the matrices ag, ,,@y, @y, 4, and ay, each of the p, being of total degree 
less than or equal to six in the elements of these matrices; and the J, are invariants 
of the form trII,y, where each II, is a matrix product formed from five or fewer 
of the factors dg, @;, Gy, Ay, As and ay. 

The matrix product ag @,@y@y 4,47 is now denoted by II, and the ij" ele- 
ments of II and y are denoted by //;; and y;; respectively. Then, from (12.11), 
(12.12), and the definition of the trace of a matrix 


ieee ic 1) 


aj Vii 
qt 
= 0] /0y;;, 
a me Pa 9 Jul OV; 5+ 
a 
But, from the definition of the f,, 0J,/@y,; is the ij element of II,. Hence 
II = a Po. tN Ree 


a 


Since the II, are all matrix products of degree less than or equal to five in the 
ap, the last relation may be written 


T= Q(5), 


which is an expression of Lemma 6. 
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Anisotropic Fluids 


J. L. ERICKSEN 


1. Introduction 

Elsewhere [/], the author formulated what may be the simplest properly 
invariant theory of anisotropic fluids. An analysis of simple shearing motion 
given there suggested that this type of theory might be applicable to some fluids 
which are commonly treated as isotropic and might describe them more accurately 
than comparably simple theories of isotropic fluids. The purpose of this paper 
is to present a revised theory which seems to us better adapted for exploring 
this possibility. This theory incorporates some features of the theory proposed 
by OSEEN [2] to describe the behavior of the so-called liquid crystals. However, 
it excludes a type of long range interaction incorporated in his theory. It seems 
to be necessary to include some such complication to obtain an adequate de- 
scription of these rather peculiar fluids. We prefer to avoid this, at least until 
the simpler theory presented here is better understand. 

Recently, GREEN & RIVLIN [3], Noi [4] and GREEN, RIVLIN & SPENCER [9] 
have formulated very general theories of materials. Existing frameworks of this 
type seem not to be sufficiently flexible to properly include the theory presented 
here. 


Cartesian tensor notation is used throughout. 


2. Conservation Laws 
Let V be any material volume with boundary S and let da; denote the vector 
element of area on S, directed outward relative to V. We assume the basic 
laws of conservation or balance of mass, linear momentum, moment of momentum 
and energy, respectively, in the forms 


(2.1) s [ eaV =0, 
4 
(2.2) 4; [ PaV = tijda;+ [ hav, 
V g 4 
4 Ss V 
(2.4) 4, [ BaV =htijkida; + [ hxav — Dada. 
Vv s Vv s 
Here 0, ~;,m;;=—m,,; and E are the volume densities of total mass, linear 


momentum, bivector moment of momentum and energy, respectively; ¢; is the 
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stress tensor, {; the body force per unit volume and gq; the heat flux. The square 
bracket means “the antisymmetric part of”, ¢.g. 
xt; =F (xf; os %;f;) . 
We envisage a fluid having at each particle a single preferred direction, 
which we represent by a vector n, of variable magnitude. One might picture 
it as a fluid composed of dumbbell molecules. Pursuing this analogy, we note 


that the total linear momentum, moment of momentum and kinetic energy, 
respectively, of a pair of mass points can be written in the form 


M<x,, M x; % 5 +Mn1;%,1, + M(4%4;%; +2; 1;) F 
Here Dou 
M=M,+M,, Mx,;=M,y,;+M22,, Mn,;=VM,M,(y;— 2), 


M, and M, being the two masses, located at points with coordinates y; and z;,, 
respectively. This serves to motivate writing 


(2.5) P= 047: 
(2.6) M5 = Pit OM Ms, 
(2.7) E= 9&+ 30(%;%;+;N)), 


where ¢ is the total less kinetic energy density, 7.e., the internal energy density. 
Equations (2.5)—(2.7) may of course be useful in describing quite different types 
of structure. 

Following OSEEN [2], we introduce one more equation of the form 


d , 
(2.8) o fear = | siav. 
V V 


He writes g; as the sum of an extrinsic part, which he suggests might be non- 
zero if there are applied magnetic fields, and an intrinsic part. In constructing 
the theory described here, we have assumed the extrinsic part vanishes. 

For sufficiently smooth functions, it is a straightforward matter to derive 
differential forms of these relations using the divergence theorem and the identity 


5 [fav = [ G+ ta, dav. 
V V 


After some simplification, one obtains 


(2.9) 0+ 0%; ;=0, 

(2. 10) o%;=t;, +h, 

(2.41) on; =, 

(2.12) Q€ = 155%; ;— ON, H;,— Q; ;, 
(2.13) +84) = by: 


If we interpret OSEEN correctly, his theory is consistent with (2.9) —(2.42), but 
not with (2.13). He assumes that ¢;;= =t;;,aSis customary in mechanics. However, 
there is some precedent for (2.13). In his theory of elastic dielectrics Touprn 
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(6, eqn. (9.10)] introduced an equation formally identical with (2.13), with 7, 
interpreted as the polarization vector, g; interpreted as an intrinsic local part 
of the electric field and ¢;; as a local part of the stress tensor. GRAD [7, § 4] dis- 
cusses conservation laws more general than those used here. In doing so, he 


introduces an equation similar to (2.13). 


Let 
(2.14) ds = 3 (Hs, + %), 4), Wi = Hu,» %4, = Fig + iy, 
(2.15) n; =n; == O55 N;- 


The vector 7; plays an important role in the theory discussed by EricKsEN [1]. 
It transforms as a vector under time dependent orthogonal coordinate trans- 
formations, whereas 7; does not. Using (2.13), (2.14) and (2.15), we obtain 


55% 4,5 = big dag + Nig Pi; 
a es 158i) @ ij 
t,54;; + 1;8;e 
= 1, ;d,5+ aa 
Using (2.16), we can rewrite the energy equation (2.12) in the form 


(2.17) oé =1t,,d;;— gi 1; — i,j: 


(2.16) 


3. Constitutive Equations 
As a general principle to be used in constructing constitutive equations, 
TRUESDELL [8, § 19] has proposed that the stress and heat flux should be assumed 
to depend on the same variables, calling this BRILLOUIN’s principle. We adopt 
an extension of this, assuming that ¢;;, g; and g; at a particle P at time ¢ are all 
functions of the variables 
(3.1) Dp PTB Uae ee, Om 


1? eeaghe Ree 
evaluated at P at time ¢, T ee absolute temperature. OSEEN makes 
essentially the same assumption as to the variables which are relevant except 
that he does not use the principle and does allow values of 1; at particles different 
from P to enter into constitutive equations. ERICKSEN [J] deals with cases 
where the stress depends only on m; and x; ; 

The relations discussed must satisfy an invariance requirement not considered 
by OSEEN [2], but now commonly used in continuum mechanics. It implies 
that these relations must be of the same form in all arbitrarily moving right- 
handed coordinate systems. This implies that (3.1) can be replaced by 


(3.2) QO; cles N;, n;, dis, de, 
as can be established by reasoning similar to that used by ERICKSEN [J]. For 


remaining implications, we restrict our attention to the case where the functions 
considered are linear in the variables n;, d;; and T.;. Thus 


=a = Al; as Aijn iy =i Ab nm Gam ae A3,T x, 
(3.3) q;= BY + Bia, + Bind, + BiT ;, 
gi = C24 Ch A, + Cid, + CLT ,;, 
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where the A’s, B’s and C’s are functions of 9, T and ,;. In addition, we assume 
that these relations are invariant under reflections through all planes containing 
n, as well as the rotations mentioned above. It is possible that this rules out 
some fluids to which the theory might otherwise apply. In general, such fluids 
would be optically active. In any event, these assumptions imply that the 4’s, 
B’s and C’s are transversely isotropic tensors with respect to the direction 1,. 
Tensors of this type have been studied by SmirH & RIVvLIN [9], who show that 
any such tensor is expressible as a linear combination of outer products formed 
from the tensors 1, and 6;;—;”,; or, equivalently, 


(3.4) N;,  Oi7- 


For the situation discussed here, the scalar coefficients in these combinations 
can be shown to reduce to functions of 


(3.5) 0, LV =n. 

As a final simplification, we assume that 

(3.6) UPR ee od PVE Nicaea ie fr a 

when we make the substitution 

(3.7) n;,>—n;, and n;>—AH,, 

v.e. that n; and —n,; are physically indistinguishable. Again, this may exclude 

some real fluids. From (3.3), (3.6), (3.7) and the results summarized in (3.4) 

and (3.5), we find by straightforward but somewhat laborious calculation that 
bj = (Hq + Ody p+ Xe Tim Me Mon + Og iy Np) 0;; sts 


(3.8) t (4 + 5 4; _ + %g Fp mn Me Mon + Oy Ny Np) N,N; 


=p Og; + Xd, Mp; + Og 1 jp My; + Oy 2; 2; + HyQN;7;, 
(3.9) 9:= Bol; + PimeT pj; 
(3.10) Gi = (Yotrdan + V24am Mx Mm + V3 Mp Mp) 1; + 74d; 0; HY 5%;, 


where the «’s, B’s and y’s are functions of the variables listed in (3.5). Now it 
turns out that (2.13) must hold as an identity to keep the number of equations 
from exceeding the number of unknowns. This implies that 


(3.41) Ya %q— HQ, = V5 = HQ — Ay - 


4. Energetics 
To complete the theory, we follow the time-honored procedure of assuming 
a local equation of state of the form 
(4.1) é=(0,%,;,S), 


where S is entropy per unit mass, to be used in conjunction with the energy 
equation (2.17). For liquid crystals, it is customary to permit « to depend on 
nm; ; as well, as is discussed by FRANK [10]. We have the invariance requirement 
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that (4.1) retain its form under all rotations, which means that it reduces to the 
form 


(4.2) é= (0, n?, S). 
We assume as usual that 
0€ 
Skee 
(4.3) il, ag = OF 


and that this may be solved for S in terms of 0, m2 and T. We may then intro- 
duce the free energy 


(4.4) g=e—TS= (0,n',T), 

such that 

(4.5) SA pe CEU byt D. 2 OU, 
00|S,n2 G0 |T,n®’ CpASioun On |tio. oar 


From (2:45), 


(4.6) N,N, =; (N; — W,;0;) =N;N,, 


since w;;—= —w,;. From (2.9), (2.12), (4.2), (4.5) and (4.6), the rate of production 
of entropy is given by 


ao Soe Op dé 
Lil BSN Ee ele 


Ses 


S = e(8— £6 — £5 2mm) 


(4.7) 


The form of the Clausius inequality most often used in irreversible thermo- 
dynamics is 


; Be a x 
(4.8) Tle S + (aslT),s] = (ti; + 8-52 O44) ds; (grt 25% ni) fy qi ;/T 20. 


More explicit relations can of course be obtained using (3.8)—(3.10). This will 
lead to some restrictions on the coefficients occurring in (3.8) —(3.10) which we 
shall not investigate here. 
5. Remarks 

Equations (2.9), (2.10), (2.14) and (4.7) yield eight equations for the eight 
unknowns 0, 7, *; and n,, it being understood that the body force and the 
relations (3.8), (3.9), (3.10) and (4.4) are specified. In general, it is consistent 
with these equations to take ;=0. The governing equations then reduce to 
those commonly used for isotropic fluids. Solutions with n;=0 are not always 
those which are most relevant. Under some conditions, solutions with »;=0 
become unstable and there are stable solutions with m;=+-0. Our investigations 
of this situation are too incomplete to be included here. 

Under the assumption that the “molecular inertia” 97; is negligible, we can 
replace (2.11) by the equation g;=0, where g; is given by (3.10). Ordinarily, 
this can be solved for #7;, which yields equations of the form 


(5.1) nN; => h;(o, ie My» Imp) - 
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Using this, we can eliminate 7; from the stress. ERICKSEN [J] introduces equations 
of this type, with @ and T omitted, on the grounds that they seem to be the 
simplest type of equations which might be used to determine ;. In this approxi- 
mation, (2.3) yields t;;=t;;, as is assumed by ERICKSEN. Generally, inertia of 
suspended macromolecules or macroscopic particles is neglected in mechanistic 
theories of suspensions. 

TRUESDELL [8] uses dimensional analysis to obtain restrictions on consti- 
tutive equations. This can yield additional restrictions in cases where one has 
a priori knowledge concerning the types of dimensional constants which can 
occur in the equations. Molecular theories might supply this. For the theory 
discussed here, one might appeal to molecular theories of suspensions of dumbbell 
molecules such as have been discussed by PRAGER [11] and GrEsEexkus [12]. 
Except perhaps in certain limiting cases, there is not much similarity between 
the two types of theories, so we shall not pursue this. We impose no restrictions 
on the types of dimensional material constants which may occur. This means 
that the form which the basic equations take for some fixed choice of units is 
not restricted by requirements of dimensional invariance. 

No unusual problems are encountered in modifying the present equations to 
obtain a simpler theory to describe isothermal flow of incompressible fluids. 
For this, we ignore (3.9) and (4.7) and replace the coefficient of 6;; in (3.8) by 
an arbitrary scalar — p. Terms proportional to d,, will of course drop out because 
of the incompressibility condition. If we also adopt the approximation leading 
to (5.1), the governing equations reduce to the form 


(5.2) bj=tj = — 2 Oj; + (Art Asda m Me Mm) Ni Nj + Ag4i 5 tq (4, M,N; +4; _ M,N)» 
(5,3) 15 = (My Hb Me Tj My Mm) My + Mg ds 52; 

(Sie PRES Pee 0) 

(5.5) te pth = 04;. 


Here the /’s and w’s are simple functions of the «’s and y’s introduced earlier. 
For fixed @ and T, they thus depend only on ®. The theory of ERICKSEN [7] 
treats m; as a unit vector, which leads to further simplification. 

To obtain properly invariant theories of anisotropic fluids, one must permit 
the stress to depend on something besides d;;, as is noted by Nott [13, § A]. 
The simplest possibility seems to be to adjoin a vector, as is done here. In some- 
what similar theories of solids, triads of vectors have been introduced by several 
writers. Information concerning this is given by ERICKSEN & TRUESDELL [14]. 
Introducing more vectors or higher order tensors would enable one to incorporate 
in theories of fluids more detailed descriptions of complex molecular structures. 
For polyatomic gases, it is conceivable that theories of this type might be useful 
in describing behavior observed in ultrasonic experiments. As can be seen from 
the work of TRUESDELL [15], the Navier-Stokes equations seem to be rather 
successful in predicting the behavior of monatomic gases. 
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The Distribution of Thermal Stress in an Infinite 
Elastic Solid Containing a Penny-Shaped Crack 


7. OLESIAK & I. N. SNEDDON 


Communicated by E. STERNBERG 


1. Introduction 


In recent years there has been great interest in the distribution of stress in 
the neighbourhood of Griffith or penny-shaped cracks in the interior of an 
elastic body, mainly because of its importance in the theory of fracture. An 
excellent review of these calculations is given by IRwrn (1958). Although, at 
the same time, there has been a revival of interest in the calculation of thermal — 
stresses (GOODIER 1958), the problem of calculating the thermal stresses in the 
neighbourhood of a crack in an elastic body does not appear to have been con- 
sidered. In this paper consideration is given to the calculation of the distribution 
of stress in the neighbourhood of a penny-shaped crack of the type considered 
previously by SNEDDON (1946). Similar calculations have been begun on the 
distribution of stress in the neighbourhood of a Griffith crack, but these will be 
reported in a separate communication. It is assumed that the stress is set up 
by the application of heat to the surfaces of the crack. 

The use of the classical equations of elastic equilibrium restricts the applica- 
tion of the analysis to materials which obey Hooke’s law fairly closely, and to 
thermal disturbances which lead only to infinitesimal strains. These disturbances 
are also required to be constant with respect to the time. By making the further 
assumption that the thermal conditions on the upper surface of the crack are 
identical with those on the lower surface, we are able to reduce the problem to 
one concerning the distribution of thermal stress in a semi-infinite elastic solid 
and therefore to make use of the methods developed in the earlier paper (SNEDDON 
1946). The solution so obtained will give an accurate picture of the distribution 
of thermal stress in an elastic body which is so large that its linear dimensions 
are very much greater than the radius of the crack. For simplicity we assume 
that the body is isotropic (though the method can easily be extended to the 
anisotropic case) and homogeneous. 

In §§2 and 3 the solution of the statical equations of thermoelasticity is 
derived for problems concerning a semi-infinite elastic solid in which there is 
symmetry about an axis normal to the boundary plane of the solid. In § 4 this 
solution is used to determine the distribution of temperature and stress in an 
infinite elastic solid containing a penny-shaped crack of radius a, across whose 
surfaces there is a prescribed flux of heat Q(7/a) which is a function of the radial 
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distance 7 from the centre of the crack. Two special cases are considered in detail 

in §5. In the first case Q(7/a) is a constant. In the second case it is assumed 

that the function Q(7/a) is such that it possesses a Fourier-Bessel expansion of 

the type >) Q,, Jo(jm7/a), where 7,,72,... are the positive zeros of the Bessel 
7 


function Jy(x). In this latter case we obtain the striking result that the normal 
component of stress is of the form Kf (r/a, z/a), where 7 and z are cylindrical 
coordinates and the function f is independent of the Fourier-Bessel coefficients 
Q,, Qo,... of the flux function; the constant K does however involve these 


coefficients. It can in fact be shown that for flux functions possessing such an 
1 


expansion that KK is proportional to [{ s Q(s) ds. A table of values of the function 
0 


f is given and its variation with 7 and z illustrated graphically in Fig. 1. This 
unusual behaviour of the normal component of stress o, is reminiscent of the 
degeneracy first noted by STERNBERG & McDowELL (1957) in their analysis of 
the thermal stresses in a semi-infinite solid whose surface is free from stress; 
they found that a, is identically zero at every point of the solid. 

Finally we consider the distribution of stress produced by prescribing the 
temperature at the surface of the crack. By calculating the equivalent flux 
function necessary to produce this prescribed temperature, we reduce the problem 
to that considered in § 4. To illustrate the method, we derive the solution for 
the case in which the surface of the crack is maintained at a constant temperature. 


2. Solution of the Equations of Thermoelasticity 


We shall consider the temperature and displacement fields in a perfectly 
elastic solid which is conducting heat. With regard to both its mechanical and 
its thermal properties the solid will be assumed to be isotropic and homogeneous. 
In the problems which we shall consider in this paper it will be assumed that 
there is symmetry about an axis, which, for convenience, we shall take to be the 
z-axis. The position of a typical point of the solid may then conveniently be 
expressed by a triplet of cylindrical coordinates (7, gy, z). For symmetrical 
deformations of the solid the displacement vector u of such a point may be taken 
to have components (u,0,w) in this coordinate system, and the only non- 
vanishing components of the stress tensor will be o,, o,, 6, and T,,. 

The equations of equilibrium of an elastic solid conducting heat may be 
written in this coordinate system in the forms (SNEDDON & BERRY 1958, p. 125) 


Oru 1 Ou u rl Vee dak ha eek oo 
pega (Sed Y OY Ay ae ») az2 ' Or az (15 2) Or’ (1) 
Cw She CIN. sh ew 0 iS *) A oo 2 
¥ 2») ( or? i y a "ra reria Oz . Oz \ oY z Y AUN OF ~ (2) 


where #@ is the deviation of the absolute temperature from the temperature of 
the solid in a state of zero stress and strain, « is the coefficient of linear expansion 
of the solid and y is its Poisson ratio. The temperature field in the solid is 
determined by Laplace’s equation 

eo 1 88 , Oo _ 


Say ir ee (3) 
in the steady state, and in the absence of thermal sources. 
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We can readily reduce this set of partial differential equations to a set of 
ordinary differential equations by introducing Hankel transforms of the in- 
dependent variables wu, w and #. We define the first order Hankel transform of 
the radial component of the displacement vector by the equation 


a(E,2) =frutr2) Aen ar, (4) 


and the zero order Hankel transforms of the z-component of the displacement 
and the temperature by the equations 


=frv(r,2) JolEr) ar, (5) 


9 (, 2) =froer, z) Jo(E7) dr. (6) 


If we multiply both sides of equation (1) by 7 J,(§7), integrate with respect 
to y from 0 to oo, and make use of well-known properties of Hankel transforms 
(SNEDDON 1951, p. 60), we find that this equation is equivalent to the second 
order ordinary differential equation 


[4 — 2v) D?— (2 — 2») 2]a —EDD= — 2(1+7) a€d, (7) 


in which D denotes the operator d/dz. Similarly, if we multiply both sides of 
equations (2) and (3) by 7 J)(€v) and integrate over all positive 7, we find that 
these partial differential equations are equivalent to the pair of ordinary dif- 
ferential equations 


[2(4 —») D?— (1 — 20) 2] 9+ EDH =2(1+7)aD9, (8) 

(D?— £2) 8=0. (9) 

If we eliminate “ and w from equations (8) and (7) in turn and make use of 
equation (9), we find that 

(D287) 420 pre Dta1b4)* wi=10. (10) 

The boundary values appropriate to any given problem can similarly be trans- 

formed to boundary values on the Hankel transforms #, @ and #, and the equa- 

tions (8)—(10) solved to give the expressions for these Hankel transforms. From 


these expressions we obtain the required solution by making use of the Hankel 
inversion theorem (SNEDDON 1951, p. 47) to obtain the equations 


)= SEE 2) Wena, 4a) 


fw (5,2) Jo(Er) dé, (11b) 


EO(E, 2) Jo(Er) dé (11¢) 
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3. Solution for a Semi-Infinite Solid 


In the case of a semi-infinite solid z=0, assumed free from disturbance at 
infinity, we are interested in solutions of the equations (1)—(3) which tend to 
zero as > oo, so that the appropriate solution of equation (9) is 


O(E, 2) = A()e"™, (12) 


where A(é), which is a function only of the parameter é, is determined from the 
boundary conditions of the temperature field. If we substitute this expression 


for & into equations (7) and (8), we obtain the pair of simultaneous ordinary 
linear differential equations 


[(4 — 2) D?—- 2(1 — ») @)a@— EDD = — 2(1+7) aFA(é) ec *, 
[244 — ») D?— (1 — 2r) 2] o+éD a = — 2(14+») aE A(d) 8, 


which are easily seen to have solutions 


UA PP er BR, ye | w= (A, bez Bye (13) 
where 
B,=B,, A,—A,+ (3 —4y) Bp =2(1+7)aE4A(é). (14) 


To find the further necessary relations among A,, B,, A, and B,, we shall 
have to make use of the boundary conditions. The shearing stress t,, is deter- 
mined in terms of the displacement vector by the equation 


Peele eee © 
A Bel Ne PRLIBE YF 


If we multiply both sides of this equation by 7 J, (€7) and integrate with respect 
to v from 0 to oo, we find that 


rt, IulEr) dr =p (Dat — EB) 
0 
= —wé[A, — B+ A,+ (B+ B,) Ezje**. 


The boundary condition 


therefore implies the relation 
Us 32 Breeden 25s (16) 


From the three equations (14) and (16) we find that 


1 + ar? = 
ee doe 7) Aa ee me 
129 een” = 
25 2(1—») as 2(1—») Sees ue) 


Now from the Duhamel-Neumann stress-strain relation (SNEDDON & BERRY 1958, 
p. 123), we know that 


1—v Ow v ou uU 4+» 
eat a (s+) aot, 
1—2v 02 1—2y \ or Y 1—2» 
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from which it follows that the Hankel transform, 
co 
G,=J 1 Jo(Er)o,47, (19) 
0 
of the normal stress o, is given in terms of “7, @ and 8 by the equation 


é,= 72 _[(1—») Dw + EH — (1+) a9]. (20) 


Substituting from equations (12) and (13) into this relation, we find that 


6,= — 74g e*{(1—») (Ag By) —9 A, H(1+9) 2 £4 At [(1—9) Ba» By] Ez}, 


and, making use of equations (17) and (18), we see that this is equivalent to the 
relation 


G,=——" [Ant (149) 2644] (1+ &2)be7*. 21) 


Inverting the expressions (12), (13) and (21) by Hankel’s inversion theorem, we 
obtain the expressions 


ster, | tes — 2»)]Ag(E) + (1+) (1+ €2) aE 4A (AF e* J, (Er) EE, 
w= f {4s) + sey [Aalé) + 1+») ae A(e)]} em JulEr) a8, 


Gf BAe JolEr) dé 


for the components of the displacement vector and the temperature. We note that, 
for this solution, when z =0, 


T= 9, as [et (1+) £7 4(E)] Jo(ré) dé 


It is convenient to write ee solution in a dimensionless form. Putting 
(1+)«A() =ap(Ea), A,(é) =(a/E)p(aé), where a is some typical length, and 
changing to a new variable of integration 7 = £a, we find that the solution we 
have obtained may be written in the form 


san |e )+ (4+) p(n)} Si(en)e*%dy, (22) 


[e,2) 


{vl + sc! [vn + 9()]} Jalon) edn, (23) 


Bete ln ~(") Jo(on)e*" dy, (24) 


where 
e=ra, C=z/a. (25) 
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This solution has the property that, when z =0, 


w =p) ) Jo(en) di, (26) 
8 Teg | nen) Jo(en) dn, (27) 
0 
a0 
Fe awa TP Solon) a, (28) 
C=O" . (29) 
etme pare < [vive ) + p(n)] Joon) ay. (30) 
Combining equations (30) and (27) we see that, for this. solution, 
o,=— Ae { nv(n) Jolon) dn — + ano, On — 0a 1) 
In general, we have ; 
a= — GH f [etn + elt +n6) e* Jolon) dn, (32) 


The determination of the unknown functions w(7), y(7) can be achieved only 
if we further specify the mechanical and thermal conditions on the surface of 
the solid. 

For instance, if the plane boundary z=0 is free from applied stress, then 
T,,=0,=0 on that plane, and we see that we must take y(7) =— p(y), so that 
the appropriate solution is given by the equations 


aor ) Js(on) eo" dy, 
0 


w=—f v(y) Jolon)e >" dy, 


0 


a wae” (7) Solon) e-*" dn. 
0 


The function g(y) is determined by the thermal conditions. For example, if, 
on z=0, it is prescribed that #=0, /(r/a), we find that w(n) =(1+7) «dy af(n) 
where /(y) is the zero-order Hankel transform of the function f(g). Hence we 
obtain the solution 


w= (1+9 ) Boa Si ) Ju(on) edn, 
w=—(1+r)ad va STty ) Jo(on) ete: 


This is the axially symmetrical solution of the problem considered by STERN- 
BERG & McDowELL (1957) and by SNEDDON & LocKETT (1960), and, as proved 
by these authors, we find that for all values of 7 and z, ¢,=1,,=0. 

16* 
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4. Prescribed Flux of Heat across the Surfaces of a Penny-Shaped Crack 


In this section we shall suppose that a penny-shaped crack, situated in an 
infinite elastic solid, is opened out by the application of heat to its flat surfaces. 
We shall assume that the flux of heat across the upper surface of the crack is 
exactly the same as that across the lower surface, and that its dependence on 7 
is known. We suppose that the crack is created in the interior of an infinite 
medium and that it is penny-shaped, occupying the circular area y<a in the 
plane z=0. If we assume that the crack is deformed by the application of an 
axially symmetrical pressure (in addition to the heat flux), then we see that, 
when z=0, 6,=—P(r/a), a prescribed function of v/a defined for y<a. As in 
a previous paper (SNEDDON 1946), we see that. we may assume that the distribution 
of stress and temperature in the neighbourhood of a crack in an infinite elastic 
solid is identical with that produced in a semi-infinite elastic medium 720, 
when its boundary surface z=0 is subjected to the boundary conditions 


oe tian O<r<a, (33) 

oo 

5, =o >a, (34) 
(0), TESA; (35) 
0,=—b(r/a), OSr<a. (36) 


From equations (33) and (34) and equation (28) above, we find, on applying 
the Hankel inversion theorem, that the appropriate form of the function (7) 
is given by the equation 


n @(n) = — a rayora® {010 (@) Jo(on) de 


which we may write in the form 


P(n) = — (1+) «a? n*9(n), (37) 
where 
q (1) aes e (2) Jo(en) de. (38) 
It follows from equation (24) that 
= sent) Jo(en) e'"dy. (39) 


Substituting from equations (36) and (39) into equation (31), we see that, when 
2=0, 0<0<1, 
_ wae (Q+)aae fr 
E~ fn) Jo(on) dy + Cree fac) Jo(on) dn = — p(Q). 
0 . 0 
Hence, if we write 


Thermal Stress around a Crack 245 


and make use of condition (35) with equation (26), we see that y(7) must satisfy 
the dual integral equations 


Savin) ) Jo(on) dn =f(e), O<e<1, 


Jol ) Jo(on) dn = 0, o>1. 


The solution of these equations is known to be (SNEDDON 1951, p. 65) 
1 Ss 
eg at * xfla)dx 
y(n) = - fsin(n s) ds Vet) (40) 
0 0 


Now, interchanging the orders of integration and making use of standard properties 
of Bessel functions, we find that 


es fau ) Jo(xn) a a= [ay [e010 Ilene 


= feos alee Jo(en) an, 


ee Be Nest Q(o)d 


Applying the rule for integration by parts to this last integral, and introducing 
the function 


0) =f e Q(e) de, (41) 


we find that this integral is equal to 


Tt 
aia Q(e) de 
220) 8 | oy) 


and, after a little manipulation, that 


1 


ee aie eyes 
v(n) =22 > f sings) ds | ae 
0 0 5 (42) 


+ (4+ naat | ao) $28". — f (9 Isler) ae}. 


We obtain the solution of the problem by substituting from equations (42) and 
(37) into equations (22), (23) and (24). 


5. Special Problems 
To illustrate the use of the results of the last section, we shall now consider 
the distribution of stress produced in the elastic solid when the flux function 
takes two special forms. In the first case we shall assume that the function 
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Q(r/a) is a constant, Qy say, and in the second that it can be represented by a 
Fourier-Bessel series. In both cases we shall assume that there is no external 
pressure applied to the surface of the crack, so that A(7/a)=0. This is not at 
all a severe limitation, since the solution for the case of a non-zero function 
p(r/a) is known for a solid with « =0 (SNEDDON 1946) and the general solution 
may be obtained by superposition. 


Case (1): Q=Qo, @ constant. 
If we put Q(0) =Q, in equation (41) we find that 


2 (9) = 2 Qo(1 — @”). (43) 
From equation (42) we find that 
y(n) = (1+) aa? Q, {nt _ 4) (44) 


and from equations (37) and (38) that 


p(n) = — (+9) 2a*Qy J. (45) 


If we substitute from equations (44) and (45) into equation (22), we obtain the 
expression 


wa Cine ee (429) \ ly? Ji ln) = 3977 cosa) = 
— (1+) 9? J, ()} J (on) e "dy 


for the radial component of the displacement vector. Similarly, from equation 
(23), we obtain the expression 


(46) 


oe Del COS 17] C tee 
== (147) «a? Q, K sa ioe ae cos} Jo(on)e*"dyn (47) 
for the other ern component of the displacement vector. 


It follows from equation (47) that the normal displacement of the surface 


of the crack is obtained by putting €=0, 0<@<1 in this integral to obtain the 
expression 


iy == (1 + ¥) «a Qo eee Jo(om) dy + 


ere Hast, f (20 1) Solem ay 
Now 
(es Tie yet yn el 0<eS1, 
6 i ae si, 
and 
gf cael nh dn = [leet Q<a, 
. u o> 
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so that, integrating both sides of the last result with respect to « from 0 to 1, 
we obtain the result 


f(a TA) 


l2n ” \ Jolen) dy = — 4 log. — + (1 9). 


Using these results and the fact that 


cosh 5 + log, @ = log, [1+ (1 — ©”) ], C2045 


we find that, when 0<0<1, 


co co 


[ SR* Iolon) ant [ (Aa ——) Jolon dn =loe. (1+ Vito] +3 (0% 


0 


so that, finally, we have 


Wy = 4 (1+) wa? Qy{2log,[1+ /(1—o?)] + (1—07)}, O<e<i1. (48) 
It will be observed that, when 9 =1, w)=0, and that, when 0 =0, 
My =4(1+4) «a? Qy(2log,2 +1). (49) 


Once we have found the integral expressions for the components of the 
displacement vector, we can readily find the corresponding expressions for the 
components of the stress tensor by differentiating under the integral sign. We 
shall illustrate the procedure involved by evaluating the normal component of 
stress o,. If we substitute from equations (44) and (45) into equation (32), we 
find that 


a, = — Hh aa Ooi (2,0), (50) 


where the function f(g, ¢) is defined by the integral 


[e,2) 


t(o,) =f +6) e*"cosy Jo(on) dy- (54) 
0 
This integral can be easily evaluated to give 
He,t) =+-cost.y +4, |o cos 3p +sin3 1, (52) 
with 
2G 
R4 = (02 +C?—1)2-+4 422, taD OP arrgaeail (53) 
It should be observed that 
fap ry 
(0,¢) eS (1-+22)2 
and that 
0, Oo 
f(@,0) = 2 
(We res a 
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The variation of f(,¢) with @ and ¢ for a set of values of @ and ¢ is shown 
in Table 1. This variation is also shown graphically in Fig. 1. From the values 
of this function we can easily calculate the normal component of stress o,, by 
means of equation (50). 


a 


Fig. 1. The variation of { (0, ¢) with e and €. (a) shows the form of f (0, ¢) as @ varies, € being kept fixed, and (b) shows the 
form of f(@, ¢) as ¢ varies, 9 being kept fixed 


Table 1. Values of the function f (eo, €) 


ele’ 0.2 0.4 0.6 0.8 1 2 3 4 i lies © 


| 


O 0) 0.0148] 0.0951 | 0.2336 0.3807 | 0.5000 0.6400) 0.5400/ 0.4429 0.3698) 0.1961 
O7 0} 0.0171 | 0.1097 | 0.2578) 0.4057 | 0.5199) 0.6410 0.5391 0.4422 0.3694 0.1960 
0.4 0| 0.0305] 0.1666) 0.3407 0.4833) 0.5782! 0.6433 | 0.5365 | 0.4404 | 0.3681 | 0.1958 
0.6 0| 0.0832] 0.3214) 0.5132) 0.6149) 0.6670 | 0.6453 | 0.5320] 0.4369) 0.3659) 0.1954 
0.8 0| 0.3763) 0.7915) 0.8113 | 0.7799) 0.7656) 0.6464 | 0.5255] 0.4322] 0.3629) 0.1949 
1 


co 11.7075) 1.2185) 1.0314] 0.9161 | 0.8399 0.6396] 0.5169) 0.4263] 0.3591 | 0.1943 
1.2 |1.5076| 1.5800 1.3166] 1.0990) 0.9558] 0.8614 | 0.6281 | 0.5062] 0.4193 | 0.3546) 0.1935 
2 |0.5774|0.5844 0.5997 | 0.6120] 0.6152 0.6092| 0.5175 0.4454] 0.3812) 0.3301 | 0.1890 
3. |0.3536]0.3547| 0.3581 0.3625) 0.3679 0.3727 | 0.3769) 0.3539) 0.3222) 0.2906| 0.1810 
4 /|0.2582]| 0.2586) 0.2598 | 0.2616) 0.2641 | 0.2666) 0.2774 | 0.2766] 0.2656| 0.2495) 0.1740 
5 |0.2041 | 0.2043 | 0.2049] 0.2058) 0.2070] 0.2084 | 0.2179) 0.2226) 0.2190] 0.2124! 0.1596 
0  |0.1005|0.1005 | 0.1006] 0.1007 | 0.1009) 0.1010 0.1024] 0.1043] 0.1063] 0.1079] 0.1061 


Case (i). Q(r/a) represented by a Fourier-Bessel Series. 


If the flux function Q(@) is represented by the Fourier-Bessel expansion 


Q(@) =2 Om JoWm@)> (54) 
where 7,72, ... are the positive zeros of the Bessel function J,(z), then, by equa- 
tion (41), ae 

2 (0) a yy —m [Js (Jn) =e Ja Gm @) | , (55) 


ma1J™ 
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with 


From equations (37), (38) and (42), we find that 


@(n) = =(A + 9) aa? Jo(n Dy ae Pe a) (57) 
2 oats (jm) f{1—cosn — 1—Jo(n) , nJo(n) 
y(n) =(1+%)aa eel : ead eae (58) 


If we make use of equation (26) and the results 


re ney cosh? (1/0) — log (1/@), Or 0 <1, 
{Buenas an =| (1/@) — log (1/@) 0 
/ 0, o> 


as ) Jol e7) )dny _ Jo(jme) 
jm? ee rancae 


we see that the normal fecren of the surface displacement is given by the 
equation 


Q,—=( +9) a.a* log [1 + Y( 1 — 0?) ry ier 1(Jm) y Se letimel | (59) 
m=1 m ™m 


m=1 


If we add equations (57) and (58), we find that 
see 4 L 2 ~ Om Sx (Im) 5 60 
v(7) + v(m) (1+%) «a ay i cere ea eae (60) 


Substituting this expression into the right-hand side of equation (31), we find 
that the normal component of stress is given by the equation 


Poe Bee) (61) 

where 
G Oe Mn (jm) 62 
=) Gee 5 (62) 


is a number which is independent of o and ¢ and the function /(o,¢) is defined 
by equation (52) above. 

If we multiply both sides of equation (54) by @ and integrate with respect 
to o from 0 to @, we find that 


7 {seas 7 © Jn lino) (63) 


so that, putting 9 =1, we find immediately that 
a 
Ca rs.0(S) ds. (64) 
0 


The striking fact proved by equation (61) is that, apart from the value of 
the multiplicative constant C, the form of the variation of o, with 7 and z is 


250 Z. OLESIAK & I. N. SNEDDON: 


independent of the values of the constants Q,, Qe, --. occurring in the Fourier- 
Bessel expansion of the flux function Q(r/a). This rather special behaviour is 
reminiscent of the degeneracy noted by STERNBERG & MCDOWELL (1957) in the 
analysis of the thermal stresses in a semi-infinite elastic solid. A similar result 
does not hold for the components of the displacement vector. For instance, 
equation (59) shows that the form of w, depends upon the values of the Fourier- 
Bessel coefficients Q,, Q2,.... We can, in fact, show the relation between w 
and Q(r/a) more clearly. If we integrate both sides of equation (64) we obtain 
the relation 


ad t 
@ 0 a 


If we substitute this expression into equation (59), we obtain the equation 


wy =(1-bo) aa Clog + V(1 — 0?) | r[Ffsomas} O sop 1 Saieg: 


where the constant C is defined by equation (64). 
For flux-functions of the type 


Q(e) =o + ce? + Co*+:- (67) 
we therefore find that 
Wy = (1+) «a?{C log [1+ (1 — o?)] + W(o)}, Oe eels (68) 
where 
el = Cy 2rt2 i 
Wile tat ake oe (69) 
and 
— 1 ~ Cy 
e ery oe (70) 


For example: 


(a) Q(0) =Qo, 4 constant. 


If the flux function is a constant then it is easily shown from equation (68) 


that the normal component of the surface displacement is given by equation 
(48) above. 


(6) Q(e) =Qo(1—@). 


For this parabolic flux function we have the expression 


Wo = 4 (1+) «a? Qof{log[1+ V(1 — o?)] + (1 — o*) —3(1— oe} (74) 


for the normal component of the surface displacement. 
(c) Q(e) = Qo(1—¢%). 


For this flux function we obtain the equation 


Wo = 4 (1 +9) xa? Qo {flog [4+ V4 — 9)] + (14 — @) -3(1—o8)}. (72) 
(4) Q(e) =Qo(1— @?)?. 
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In this case we obtain the expression 


=7(1+) «a? Qo{§ log [1+ y(1 — o?)] + (40?) —4(1— 0) + $(1 —0%)} (73) 


for the normal component of the surface displacement. 


The variation of the normal component of the surface displacement for these 
four forms of the flux function is shown in Fig. 2. In this diagram, the length L 


is defined by the equation 
y q is \ x a? Qo. (74) 
0 i Q 
ib 
W, ay, 


Fig. 2. The normal component of the surface displacement in the four cases considered in § 5. In each case the upper 
curve shows the flux function and the lower curve the displacement component 


6. Prescribed Temperature at the Surfaces of a Penny-Shaped Crack 
In this section we shall discuss briefly the case in which a penny-shaped crack 
is opened out by the application of a prescribed temperature to its flat surfaces. 
In this case the boundary conditions (35) and (36) remain unaltered but equa- 
tions (33) and (34) are replaced by the pair of equations 


8=—S gla), OSrSa; i) 
Comme 
ces oe (76) 


In equation (75) % is a constant and g(r/a) is a prescribed function of 7/a. 
It is easily shown that the solution of Laplace’s equation (3) satisfying these 
thermal boundary conditions is 


jas af g(n) Jo(on) e-" ay (77) 
with 


= [na Jo(en) e” dy. (78) 


0 
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For this solution, therefore, 
Q(e)=JS a(n) Jolen)dn, 0<e<1. (79) 
0 


As a result of the boundary conditions (75) and (76) we see that the function 
q(yj) must be taken to be the solution of the dual integral equations 


A, 
fae ) Jo(on) dn =—*g(e), O<e<t, (80) 


Jnatn ) Jo(on) dn =0, eA. (81) 


The solution of these equations is known to be 


1 


con [1h 2? Bn | ihe fel) sin (qu) du (82) 


0 


(SNEDDON 1951, p. 65). If we substitute from equation (82) into equation (79) 
we obtain the expression for Q(o0), the flux of heat across the surfaces of the 
crack. Once we have found this function, we can proceed using the method of 
§ 4 above. The following results are useful in applying this method. 


From equations (41) and (79) we find that 


=S a(n ) (Sim) — e Jule) ay, (83) 
and hence that 


mi a(n) Jun) dn. (84) 


Similarly, from equations (64) and (79), we find that 


C=JSa(n) Jun) dn, (85) 
and, from equation (79) that 


o=f4 Hf 20t fae a [Jolen) — Jo(n)] dy. (86) 


For instance, if we substitute from equations (85) and (86) into equation (68) 
we obtain the appropriate expression for the normal component of the displace- 
ment of the surface of the crack, and if we substitute from equation (85) into 


equation (61) we get the expression for the normal component of stress at any 
point in the solid. 


To illustrate the procedure, we consider the simple case given by the equations 


0=—%, O57 SG, @=O05s 
oo 


eee f= Oh. Pe 0) 
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In this case g(v/a) =1 and it is easily found that 


g(q) = 200 SE 
so that 
Q(e) = 2% f sing Jo(gn) dy. 
From this expression it is readily oe that Q(e) =0 if o>1, and that 
Q(e) = 2% (1—6%)-#, ~ O<@ <1. 


UA 


By the definition of 2(0), we find that 


28 re ON 
2(0) = er (eee) 
Equation (42) then yields the expression 
ONG 0 == ~ 
Ue “a —e fva-e) Julen ade} (87) 
0 
for the function y(y). Evaluating the integral we find that, for this problem, 
_ (1+7) «Pa 4 sing 
v(n) = (2 — cosy — = i (88) 
Similarly, equation (37) gives the equation 
me 2(1+r) aa sinn 8 
y(n) ae en (89) 


The solution of the problem is then obtained by substituting from equations 
(88) and (89) into equations (22), (23) and (24). 


For this problem 
t 


fsQ(jas=2% (1 —(1—#)8], 


so that 
_ 2% 
” 32a 
and 


Le 205i 


we 37a 


Mor 0%) — leg [ty 1 — 0°). 


Substituting these expressions into equation (68), we find that, in this case, the 
normal component of the displacement of the crack surface is given by the equation 


wy = Ly [(1 — 0”)? + $1 — 0?)] (90) 


where 
L.= renee (91) 
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In a similar way we obtain the expression 


_ 2(1+r) «da 


for the normal component of stress at any point in the solid. 
The variation of wy with @ is illustrated in Fig. 3; that of o, can be inferred 
from Fig. 1. W, 
ade] 


Lo 


0 7 


Fig. 3. The normal component of the surface displacement when the surface of the crack is maintained at a temperature 
0) below the temperature of the reference state 
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Null Solutions of Partial Differential Equations 


LARS HORMANDER 


1. Introduction 
Let P(D)=P(—10/0x, ..., —10/0x") be a partial differential operator with 
constant coefficients. A solution u €&(R”)* of the equation 


(1.1) P(D)u=0 


which vanishes outside the half space H={x; (x, N)>0} will be called a null 
solution of P(D) with respect to H. It was proved by HORMANDER [1], Theorem 3.2 
that there is a null solution +0 with respect to H if and only if the bounding 
hyperplane <x, N)=0 1s characteristic, that is, 


(1.2) P(N) =0 


where P° is the pricipal part of the polynomial P. Our aim here is to investigate 
when all solutions can be approximated by null solutions. 


Theorem 1. Let P(C)=A (6) B(C) where the plane <x, N>=0 1s characteristic 
with respect to every irreducible factor of P, but non-characteristic with respect to Py. 
Let H be the half space {x; <x, N)>0}. Then the closure in &(H) of the null solutions 
of P(D) with respect to H 1s the set of all solutions in &(H) of the equation P, (D)u=0. 


In proving this theorem we may assume that P,=1. For if 
P(D) u = P,(D) (B,(D) u) =0 


and u vanishes outside H, the function A (D) uw is a null solution for P,(D), 
hence identically zero, since the plane «x, N)>=0 is non-characteristic for P,(D). 
Hence PR (D)u=0, which shows that the null solutions of P(D) and of A (D) 
are the same. Thus we may assume in the proof that the plane <x, N)=0 1s 
characteristic with respect to every irreducible factor of P(C). 


A weaker result than Theorem 1 is the following 

Theorem 2. Let P fulfill the same assumptions as in Theorem’. Then the 
closure in &(R") of the set of functions u which are null solutions with respect to 
the half space {x;<x,N >> c}, for some c, consists of all solutions in &(R") of the 
equation P, (D)u=0. 


* By &(Q), where Q is an open set, we denote the set of infinitely differentiable 
functions in 2 with the Schwartz topology of uniform convergence on every compact 
set of the derivatives of all orders, including 0. For the definition of &(Q) and 9(2Q) 
we refer to SCHWARTZ [45]. 
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2. Preliminary lemmas 
Our first lemma is of an algebraic nature: 


Lemma 1. Let Q(C) be an irreducible polynomial such that Q(sN-+&) 1s not 
independent of s for indeterminate s and &. Then the polynonual 


Q(sN +t& + &5) 


is an irreducible polynomial of s and t except when (E, &)) belongs to a certain al- 
gebraic set. Hence there is an algebraic set A such that when €¢A the polynomaal 
is irreducible except when &, belongs to an algebraic set, depending on &. 

Here we mean by an algebraic set a set defined by algebraic equations and 
different from the whole space. | 

The author has been informed by Professor P. SAMUEL that this result is 
well known though not explicitly stated in the literature. A proof may be based 
on pp. 211—218 in LaNc [3] but will not be given here. A somewhat weaker 
lemma was used in a similar context by JOHN [3, pp. 397—398]. 


Lemma 2. Let I'(C) be an an entire analytic function and P(¢) a polynomial. 
Suppose that there is a vector E such that P®(&) ==0 and 


F(tS + &o)/P(LE + So) 


is an analytic function of t for a dense set of complex &. Then F(¢)/P(C) ts an 
entire analytic function. 


Proof. Since the coefficient in P(t€+&)) of the highest power of ¢ is inde- 
pendent of &) in view of the fact that P°(€)=+-0, the zeros ¢ of the equation 
P(t§+&,)=0 vary continuously with €). For all &) in a dense set the zeros of 
P(t&+&)) are by hypothesis zeros of F(t§+&) with at least the same multi- 
plicity. Hence it immediately follows that this is true for all &). 

We can now define G(& >) = lim F(t&+ &)/P(t&+&); the limit exists because 


the quotient is an analytic function of ¢, and G(&)=F(&,)/P(&) if P(&) +0. 
Choose 7 so that 


P(tE+&)+0, |t)=r. 


Now we have by Caucuy’s integral formula 


PCR i 


lél=7 


For ¢ in a neighbourhood of & we have P(t&+¢) +40 in the range of integration. 
Hence the integral is there an analytic function of ¢, and G(¢) is thus an entire 
function. 


We finally give a lemma on Fourier transforms, which is also well known. 
The Fourier transform of an integrable function is here defined by 


@(C) =f eS*® p(x) dx. 


Lemma 3. Let 0>1 and w be a function with compact support such that with 
a constant C 


(2.1) | D* p(x) S Cal (af? 


C(t) San) ae [ae ees 
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Then the Fourier transform has the estimate 
(2.2) [P(E + ¢n)| S Cexp(A In| — BIEIM), 


where C,, A, B are constants > 0. 


Proof. From (2.1) it follows with a constant C, that 


p>) EFS O| SCH eetlml, Cae Lin, 
which gives 

, “a * C, k , 
23) PE +i] SCs) (R1)% eA lal, 


In order to make the right-hand side of (2.3)’ as small as possible we choose k 
as the largest integer <(|&|/C,)/°. Using StrrL1Nc’s formula 


bl =(Z) yank (1+ o(1)) 
we then obtain (2.2). 


3. Proof of Theorem 1 

We have to prove that every distribution TC&é'(H) which is orthogonal to 
every null solution with respect to H 1s in fact orthogonal to every solution of (1.1) 
in &(H). We first reduce our study to such distributions T as are very smooth 
functions; this is essential for the convergence of the integrals occurring below. 

Thus let # be a function in Y(—H) for which (2.4) holds and f{ @dx=1. 
Such a function exists for every 6>1. Indeed, it is well known that the function 
which equals 0 for x,;<0 and exp(—.«}~°) for x,>0 satisfies (2.1). Starting 
from this function, we easily construct a function # with the required properties, 
since the set of functions satisfying (2.1) for some C is a ring (cf. SCHWARTZ 
[5, p. 22]). Now form 7;=T* #0, where 0,(x) =e" O(4/e), €>0. When e>0 
we have T,->T in the topology of &’(H); and 7;€ Z(H) when ¢ is small enough. 
Moreover, T, satisfies the condition (2.1) for every e, hence is satisfies (2.2). 

Now it is easy to see that if T is orthogonal to all null solutions u with respect 
to H, then so is 7,. For 


vy 


T,(u) =T,* i (0) = T*9,#14(0) =TW.* 1), 


é€ 


where “(x)=u(— x). It is clear that 0% u is also a null solution with respect 
to H, hence T,(u)=0. If we can prove that 7} is orthogonal to all solutions 
u €6&(H) of (4.4), it will follow when ¢—0 that 

T (#) =tim F (uv) —=0. 


e—>0 
Hence it is sufficient to show for some 6>1 that 1f a function p€ D(H) satisfying 
(2.1) is orthogonal to all null solutions, then y is orthogonal to all solutions. 


We now recall the construction of null solutions given by HORMANDER [1], 
giving it a slightly more general form. Let 


P(E) = Qi (0)"--. (0) 
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be the decomposition of P in irreducible factors, Q,,..., Q, being different. We 
know by hypothesis (see Section 1) that 


(3.1) OF (N) 0, ieee erent 2 
Let € and &, be fixed vectors such that € is real and 
(3.2) OSE 0,2 It ieee 
It is then easy to see that the equation 

33) Q,(sN +4&+&) =0 


has a solution ¢(s) with a Puiseux series expansion 
CO 
— k/p _s lb 
t=s x Cc, Ss 


convergent when |s|=M. (Cf. VAN DER WAERDEN [6, pp. 50—52].) Further- 
more, which is the important point, we have k<# so that 


(2s) FO (iste ys iS => 00, 


Let @ be a number such that 
1—1/p<o<1, 

and set with e>0 
(3.4) u(x) = f eH sN+#(s) E+6) stl? e—elis)® dg. 

Be 
Here we define (7s)* so that it is real and positive when s is on the negative 
imaginary axis, and we choose a fixed branch of s’/? in the lower half plane. 
The contour C~ is composed by the parts of the real axis where |s|=M and the 
half of the circle |s|=M in the lower half plane. The curve is described from 
—ooto + oo. Finally y is an arbitrary integer. 

Arguing precisely as on p. 217 in HORMANDER [J] we find that wu is a null 
solution of Q;(D) with respect to H. In fact, w is a solution because of (3.3), 
and a null solution because in (3.4) the integration contour can be deformed 
to the line Im s=rt if r< —_M; when t-+— oo an easy estimation shows that u 
vanishes when <x, N)<0. The details of the proof need not be repeated here. 
Since Q;(D) is a factor of P(D) of order yw, it follows from LErpniz’ formula 
that R(x) w(x) is a null solution of P(D) if R(x) is a polynomial of degree <U;. 

Now let »¢€QY(H) be orthogonal to all null solutions with respect to A. 
Differentiating the equation 


P(E) =f eb w(x) dx, 


R(D) 9) =f eb R(x) p(x) dx 


if K is a polynomial. Thus if we multiply (3.4) by R(x) p(x) and integrate with 
respect to x, if the degree of R is <y;, so that R(x) w(x) is a null solution, we get 


(3.5) (R(D) @) (sN + t(s) € + &) see's)’ ds = 0. 


we obtain 


The change of the order of integrations is legitimate because we have an absolutely 
convergent double integral. 
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We want to let ¢0 in (3.5). In order to be able to do so we have only 
to make sure that the integral is absolutely convergent for e=0. For the factor 
exp (— e(7s)?2) is <1 in absolute value on C~. Now, as pointed out above, we 
may assume that R(D)@ satisfies an inequality of the form (2.2). Hence with 
suitable constants C,, A, B 


|(R (D) $) (sN + t(s) & + &)| 
<= C,exp (A|N Ims + £ Im ¢(s) + Im &)| — B|N Res + €Re f(s) SRE Sole ie 


Since |¢(s)|=0(|s|'—/*), it will follow that the integral (3.5) with e=0 is abso- 
lutely convergent if 


1 1 seers 
lca petasgie that is, if d<p/(p —1). 


Let m be the degree of P and take 6 so that 
1<d<m|(m—1). 


Since Pm we then have 6</(p—1) so that the integral (3.5) is absolutely 
convergent when e=0. Letting ¢«—0 we thus get, writing 


(R(D) @) (sN + t(s) & + &) = F(s*?) 
so that / is a univalent analytic function, 


(3.6) [aE stPy Hid s= 0: 
es 


On the other hand, if in (3.4) we integrate over the curve C* which is sym- 
metric to C~ with respect to the real axis and is described from + co to — oo 
and replace the convergence factor e~*/)* by the factor e~*%#*, we will get 
a null solution u with respect to —H. Hence u vanishes in the support of g, 
and is thus orthogonal to kg. Repeating the same argument as before, we obtain 


(3.7) | FC) s#? ds = 01. 
Gk 


The equations (3.7) and (3.6) are valid for all the different determinations 
of s/?. Adding the equalities (3.6) and (3.7) for all the determinations, we get 


(3.8) f F(s"?) sds =0 


where the integral is taken # times round the circle |s|=M. In fact, it is obvious 
that the integrals along the real axis cancel each other. Setting s=2?, we get 


(3.9) litte) 27a de = 0 
|z|=' 
for all integers v. Hence all coefficients in the Laurent series for F(z) are zero, 
so that F(z) =0 identically. 
Before proceeding we have to choose & ina suitable way; up to now our only 
conditions on € have been the conditions (3.2). We now also demand that for 
those Q;, 7=1,.-.,7, which are not constant in the direction N, the vector § 


? 


shall have the property given by Lemma 1. These requirements can be fulfilled 


i? 
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at the same time as (3.2) since each condition means only that & shall not belong 
to a certain algebraic set. Thus we now know that the polynomials Q; (s N + t€+&)), 
{=e push? pate utneaueible if £)¢.A, where A is also an algebraic set. 

When & ¢A it follows if 7=1,...,7’, that (R (D) ?) (sN+t&+£,)=0 if the 
degree of R is <j, and (s, ¢) is on the Riemann surface of the algebraic function 
defined by 

Q;(sN +t& + &) =0. 


In fact, this surface is connected because of the irreducibility, and we have proved 
previously that the assertion is true in a neighbourhood of infinity in a certain 
sheet. Hence it follows in general by analytic continuation. Hence p(s N +t&+ &) 
will as a function of ¢ have a u;- fold zero wherever Q;(sN+t&+& ) has a zero, 
and passing to the limit from such s for which the zeros are simple, we also find 
that a zero of Q; of multiplicity & will be a zero of @ of multiplicity & wu; at least. 
In particular, for s=0 we find that 


G (EE + §)/O; (2 + Eo)" 


is an analytic function of ¢ when &,¢ A. 

If j=r'+41,...,7 so that Q; is constant in the direction N, it is even more 
easy to get the same result. In fact, if Q;(t+&), 7>7' has a zero ¢) of order mu, 
then 

1 (1) = Cx, Eh ef mnttt 


is a solution of the equation Q;(D)“uw=0 if k<wu;u. We can find a function 
XE Co(0, 09) so that y(<x, N>)=1 in a neighbourhood of the support of 9. 
Since Q; is constant in the direction N we have 


Q;(D)i (yu) = 49; (D)u=0. 


Hence yw is a null solution, so that 


Soudx=fop(yu)dx=0. 


It follows that ~(¢&+-&)) has a zero of order w, wu at least for =t). This proves 
that s 
PEE + So)/Q; (FE + Fo) 


is an analytic function of f. 


Now let B be the set of all & such that at least two of the polynomials 
Q; (t+ 0), 7=1,...,7, considered as polynomials in ¢ have a common zero. 
B is an algebraic set since it is the set where the resultant of Q; and Q, vanishes 


for some pair 7==k and none of these resultants is =0. When &)¢(AUB) we 
now get that 


PLE + La)/PUCE + £0) =O (68 + 60)/TT Q,(tE + So)" 


is analytic. In virtue of Lemma 2 it follows that $(¢)/P(C) is an entire function. 
Hence according to MALGRANGE [4, Chap. I, Prop. 2], we have @ (C)/P(C)=(C) 
where yp €&’(H) (in fact, y€Q(H)). From the equation @ (C) = P(C) p(C) it follows 


that 
p=P(—D)y. 
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Hence if w€ &(H) and P(D)u=0, we get 


fu pdx =(P(—D)y) (w) =y(P(D) «) =0, 


which completes the proof of Theorem 1. 


4. Proof of Theorem 2 
Again, what we have to prove is that if @ is orthogonal to all null solutions, 


then is orthogonal to all solutions of (1.1). But if we take the half space H 
so that m has compact support in H, it follows from what we have already proved 
that p= P(—D)y where p has compact support in H. Hence @ is orthogonal 
to all w€ &(R") with P(D)u=0, which proves the statement. 
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The Formulation of Constitutive Equations 
in Continuum Physics. II 


R. S. RIVLIN 


1. Introduction 


In previous papers by GREEN & Riviin [/, 2] and an earlier report [3] 
constitutive equations have been developed describing the mechanical propertie 
of materials possessing memory. A similar topic has been discussed by Nott [4, 3]. 
In [J, 2] it is assumed that the stress at a fixed particle of the material is a func- 
tional of the displacement gradients at the instant of measurement and at times 
preceding the instant of measurement. In another paper, by PrpKIN & RIVLIN [6], 
constitutive equations have been developed governing certain types of physical 
properties of materials which may be described by the relation between a tensor 
of order yw on the one hand and the displacement gradients in the material and 
y vectors on the other, all the quantities being measured simultaneously. In the 
present paper, these two problems are combined. 

We are concerned with physical properties of materials which can be described 
in a rectangular Cartesian coordinate system by the relation between the com- 
ponents at some instant ¢ of a tensor u of order wz on the one hand and the com- 
ponents of » vectors and the displacement gradients on the other, at all times 
from some fixed initial time up to and including the time ¢. The components 
of wu are therefore functionals of the » vectors and of the displacement gradients 
over a time interval preceding ¢. 

The limitations imposed on the constitutive equation by the assumption that 
it is unaltered by a simultaneous rotation of the reference system and of the 
physical system are derived. The manner in which we may derive the further 
restrictions due to the symmetry of the material is discussed and illustrated in 
the particular case when the material is isotropic and possesses a center of 
symmetry. 

Finally, in §§6 and 7, we discuss the manner in which these constitutive 
equations may be considerably simplified if we restrict ourselves to certain 
classes of deformations in which the displacement vector and » independent 
vectors in the constitutive equation each vary in a prescribed fashion with time. 
§ 6 deals with the case in which the body is undeformed and § 7 with the case 
in which the displacement gradients are small compared with unity. The 
approach taken is similar to that adopted in a previous paper [7] in which is 
discussed the constitutive equation appropriate to purely mechanical problems 
in which the stress at any instant of time depends on the displacement gradients 
at and preceding that instant, the displacement gradients are small compared 
with unity and the displacement vector varies in a prescribed manner with time. 


. 


n 
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2. Restrictions on the constitutive equation due to rotation in implicit form 


We consider a body to undergo deformation such that a generic particle of 
the body which at time t=0 is at X; in the rectangular Cartesian coordinate 
system x moves to x; in the same system at time t. Then, the deformation of 
the body is described by the dependence of x; on X, and t, thus: 


x, = %;(X,,T). (2.1) 


We assume that some physical property of the material is described by the 
dependence of the components 4; ;,__;,(f) of a tensor of order at time ¢ on the 
values of the deformation gradients 0x,(r)/0x, and of the components in the 
system x of y vectors v(t), (x =1, 2,...,”), at all times in the interval 0 to t. 
This may be expressed mathematically by 


t 
1 ae ieee LIES oe) , (2.2) 


where fj ;,..;, denotes functional dependence on the arguments. In (2.2), we 
employ the notation 
en eee aye (2.3) 


The assumption will be made later that fj, ;,_.;, is a continuous functional of its 
arguments in the interval 0 to t and is particularly dependent on the values 
of the arguments at t =1. 

We assume also that the related physical variables are such that the con- 
stitutive equation (2.2) is unaltered by a rotation of the physical system as a 
whole, if it is expressed in a coordinate system which is subjected to the same 
rotation. The restrictions which are thereby imposed on the constitutive equa- 
tion (2.2) may be obtained in the following manner. 

Let us assume that in addition to the deformation described by (2.1), the 
body is subjected, subsequently to time t =0, to a time-dependent rigid rotation 
in which a point at x;(t) in the system x is moved to %,(t) in the same system, 


U 


where a 
%,(T) —— a; ;(T) Xj (t) c (2.4) 
The quantities a;;(t) satisfy, of course, the orthogonality conditions 
4; ;(T) 4;,(T) = 4; ;(T) a; (T) = Ojn, |a;;(z)| =1. (2.5) 
We then have directly from (2.4) 
O%5(t) 64, (7) m 
aX, - = Ay, ( ) OX, a (2.6) 
We consider also that each of the vector fields ae (a) (o==4) 250... v7) 1s, subjected 
to the same rotation. Then their components DONT) in the system x are given by 
Be Garros 0) (2.7) 


Suppose that if the tensor field w; ;,;,,(7) is subjected to the same rotation, it 
has components #;,;,...;,,(t) in the system x. Then, 


Ui, ig... Mii, oc ] uy, I ] (2.8) 


where the notation ba 
a; ; — a; ;(t), uj, ty... OU = Uj, i,...in(t) (2.9) 


is used. 
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Let x’ be the coordinate system into which the coordinate system x is rotated 
by the rotation described by (2.4). From our assumption that the constitutive 
equation is unaltered by a simultaneous rotation of the reference system and of 
the physical system as a whole, we see that in the coordinate system x’ the 
constitutive equation is (2.2). In the coordinate system x the constitutive equa- 
tion may, of course, be written as 


Ui iy...i =F; j, 


t 
: iol 2F (2/2Xy, BY (0)]- (2.10) 


From (2.2), (2.8) and (2.10), we obtain 
t 
Fig-in| 2% (EXg, BHC 2) 


t= 


t 
| a Wis Bigjgs** Fivin ti, fei OX) oX4, u(r) (2.41) 
This gives in implicit form the restriction imposed on the functionals Fj, ;,. i, 
by invariance of the constitutive equation under simultaneous rotation of the 
physical system and the reference system. In the next section, we shall obtain 


these restrictions in explicit form. 


3. Restrictions on the constitutive equation due to rotation in explicit form 

We replace the functional dependence of F;, ;,.._;, in (2.2), and hence in (2.11), 
on its arguments 0x, (7)/0X, and via), by polynomial dependence on the values 
of 6x,(t)/OX, and v(t) at N+-14 discrete instants of time tp (P =O, 1, 2, ..., N) 
in the interval 0 to ¢, where t,=t. Thus, 

Ui i,...ip Leh [ox, (tp)/OX,, oS (tp) | , (3.1) 
where fj, ;,...i, are polynomial functions of the indicated arguments which satisfy 
the relation 
Ve sp [Ox,(tp)/OX,, o. (eel = Gi, 5j, Gigjg- ++ Vin ju hi. jes-- ips [Ox,(tp)/OX,, vu!) (tp) ]. (3.2) 
This relation must be valid for all a, ,(t) satisfying (2.5). 


Let us suppose first that a;;(t) =0, ; at all instants of time except in a small 
time interval about tg =+-¢. Then (3.2) becomes 


ee OX p(T) AX (Ty) 0% (TQ) 0x p(T) 
Vy y..-4 lb ta 5) 1 iy aeRel:y ye Deh) 
Blane aX, OX, Pha 
Oe (ae), DBCAE Saf oe (TG) Aan; 0) 
a OX p(T)  Oxp(T,) 0X» (To) OX» (TN) 3:3) 
tte tp) ay , pee MIX BT BCS 
aX, aX, aX, Sa. ONe. 
UO) Do (,); toe vi)(t9), ata: vi (ty) ; 


This implies that Tel gok ity is a polynomial scalar invariant of the y+ 3 vectors 
OXp(To)/OXy, OxXp(T)/OXy, Oxy (T)/OXg and v (to) (x=1, 2,...,”). It follows 
that Vigig. vi must depend on OX (To)/OX, and vs") (Zo) through the elements of 
an integrity basis for the » vectors. Such an integrity basis is given by 


Gas (To), VAG als wee) o}(c9) 5) 


0x;(T) Ox; (tr, 
[G(t9) ]2, Ci ik ax, aS gue (to), 
q 


045 (TQ) (a) 
ao em Y(t) W(t), ej V(t) vP(zg) Ca 
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where the notation 


OX p(T % 
Gra 0) = oe) SR) GQ) =(G.(0)] and Y% er) = 80) olor), (3.5 
is used and e;;,=1 or —1, accordingly as 7,7, k is an even or odd permutation 
of 1, 2,3 and e;;,—0 otherwise. 
In a manner similar to that employed in a somewhat similar problem by 
PIPKIN & RIVLIN [6], it can be shown that any polynomial in the quantities (3.4) 
may be expressed as a polynomial in 


G,s(t), [G(t)]-# and V'"'(t9). (3.6) 


Repeating this argument for each of the instants of time 1, T),..., Ty, We oe 
at eq. Ga musts be pxmeeable as a polynomial fh; ;, an ae in Ax pK, oy , 
Gy 4 (Tp); (G( (Celle tr? (ee 12 at I) l= 1) 2a t-ad) thes: 


On. 
lita, nies ve); Gy 4( Tp), ), {G (tp)}- , Vie) ) 


= Hive inl BX, 


(3.7) 


ty 


where the notation x,=~,(t) and v ae ean) is used. Introducing (3.7) into (3.2) 
and taking a,,;(t) =0;,; for all instants of time except in a small time interval 
about t =¢, we obtain 

OXp (a) 


Ox. 
p 

ee, b 

Ox pat 


Die 


(a) 


(3.8) 


ee of h,, i i, on both sides of (3.8) on the remaining variables G,, (Tp), 
[G(tp)|~# and ye (tp) being understood. It has been shown in the previous 
paper [6] that a Sete Con of the form (3.8) implies that /,,;,.__;, must be expres- 
sible in terms of ¢x,/0X, and Ney by 


|= = Gi, 5, Vinings ++ Gi ju My jy. 


OX, OX}, OXin ( 
Nii,...in = OX; OX), . Tee, Pj, jo... Gage Cus ale (3.9) 


where 9j, ;,,..j,,.18 @ polynomial in the indicated arguments and in the arguments 
Gy (tr), [G(tp)]~4 and Vj (rp), (P =1, 2,..., N). From this result together with 
(3.7) and (3.1), we obtain 


Ox; O%;, ax 4 7) 
Ui, i9...ip 0x, OG cea aX, Pas joes iu [Gyq(tp), {G (tp)} a) V5" (tp) |, (3.10) 


where P =0,1,2,...,.N. Since we may express {G(tp)}~? as a polynomial in 
G,,(tp) to any desired degree of approximation, we may without serious loss of 
generality omit the arguments {G(tp)}~? in (3.10). Now, allowing the number 
of instants N to increase indefinitely, we obtain 

OX;, OX; Ox; } 


mee eat UBUD ee (a) 
Ui in... ip aX;, aX;, ae OX}, Di Gr Coal), Vo (T) ? (3:1) 


where ®, ;,.__;,, iS a functional of its arguments. 


4. Restrictions on the constitutive equation due to material symmetry 

We shall now consider the further restrictions imposed on the constitutive 
equation (3.11) by any symmetry which the material may possess in its initial 
state, in which it is undeformed and of) =0. The constitutive equation must be 
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form-invariant under the transformations from one rectangular Cartesian co- 
ordinate system to every equivalent rectangular Cartesian coordinate system. 
If x and % are equivalent coordinate systems related by 


x(t) = S54 x; (T) , (4.1) 
where Sij Sip = Sji Shi = Oj ke (4.2) 
and the constitutive equation takes the form (3.11) in the system x, then it 
takes the form 
= O%;, 0%; Xin @ ras yi) 
gu, See feta ee a wpa Gs 4.3 
Ui ig...ip aX;, axy, OX jy Nidow- 7h bah )s Vp ( ) ( ) 


in the system ¥, where %; is given by (4.1), 


X;=$;;X;, Gp q(t) = Spy Sq5G,5(T), (4.4) 
Ae = Spq ee and Ui, to... tp es Sis; Sinise ee Sin eu Uj, Jo-+- Tp 
From (3.11), (4.3) and (4.4), we obtain 
OX;, OX; OX jy | : r(x) 
So Sia ke Shay z ae ENG: G Gy), VAG 
Sia; ints Siu ju axe Oxy. OX py Ry ky... Ru veal ), ( ) 4 5) 
O%;, 0%), Fig = t 


OX}, aX), OX}, Tifa Th 
This restriction on the form of the functional @;,;,_;,, must be valid for the 
transformations relating the coordinate system x with each of the equivalent 
coordinate systems. We can obtain this restriction in a more explicit form in 
the following way. 

Let wp”) and py) (6 =1, 2,..., uw) be the components in the systems x and % 
respectively of w vectors. Then 


Pr) = 85.9 (4.6) 
pe ee pp?) ple) i i 
Multiplying (4.5) throughout by »;)’ y;?)...p?, employing (4.6) and the notation 
B J Ox 7 aS Ox = 
yA) — ylP) Mee PFA) — a0) x (4.7) 
we obtain 
t 
FAN A A Ds. n,.. in| Gog (t), VEO) 
=0 
= (Me Aap Yio. UE Stet [G,,(t), clic (4.8) 
=F (say). 
With this definition of F, we note, from (4.7) and (3.411), that 
oh On 
Wiig... ip dy dy... dy) : (4.9) 
Vy ty ue 


We now assume that, in our initial equation (2.2), the functional F& iinds 
a continuous functional of OX» (t)/OX, and v(t) («=1, 2,..., y), over the compact 
aggregate of these functions which are continuous in the range OStS/, and is 
particularly dependent on the values of these functions at the instant ¢. It 


follows that the functional DN ai is a continuous functional of G,,(t) and 
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V(x ), over the compact aggregate of these functions which are continuous in 
the range 0S tSzt, and is, in general, particularly dependent on the values G, , 
and ae of these functions at the instant ¢. We assume sufficient regularity = 
that ine , may be expressed as a polynomial in G,, and V{, the coefficients 
of which are continuous functionals of Gyq(t) and "Ye Tt), over the compact 
aggregate of these functions which are Eeauiaiene in the range OS tS. 


We shall denote the Fourier half- as cosine coefficients of G,,(t) by yl’? 
ea Oul 229 se )hand Whose: of "VA"'(7) by “Oa P=05.1 254. 5) where 


2 IP 
0 =2 f Gygl0) cos Sites lla Aih 
t 


’ (4.10) 


and ; 


Then each continuous functional in the expression for ®; ;,.__;,, may be approxi- 
mated by a polynomial 
(N)/,,(0) (1) (N) (0),(a) (1) 
Pe (year Voge -291 Yoqn Op» D> 


which tends uniformly to the functional as N— oo. It follows that ®;,;,_;, may 


be expressed to any desired approximation as a polynomial P*,, ain 
Pp y PP poly ps 
(0) ,,() (N) (0). (a) (1),,.(a) (N) (a) 
ca poe Peden pe, Ob A ers HE 


Sa ee.) Ne and “)@\” in terms of G, q(t t) and ye (t) in the same way as Na) 
a) 


and | oe are ae in terms of G,,(t) and V{"(t), we have 
malta 
Yoo = Soest. 


(PEO) — 5, Pen, 


and (4.11) 


We note also from (4.1), (4.4), (4.6) and (4.7) that 
Fe) — 5, we, (4.12) 


Replacing ®,;,..;, in (4.8) by Ba iz» it follows that F is a polynomial scalar 
invariant, under the group of transformations characterising the symmetry of 
the material, of the w+v(N-+ 2) vectors TAG (Cs Mee UA Vis (oped Di oe V) 
and 228 Marae 2,..-,¥7; P=0,1, 2,..., N) and the N-+-2 symmetric tensors 
G,q and Yea (P==0)1,2,-.5, NV) It iollows thate® mustbe: expressible asea 
polynonial in the aimichts: oi an integrity basis for ioe vectors and tensors 
and we note that F is multilinear in the vectors 1 (i= 42,5, 4), sowthat 
we need not consider those elements of the Atco basis which are non-linear 
in any one of the vectors Y”. 
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5. Isotropic materials possessing a center of symmetry 

Following the argument of the previous section, we see that if, prior to time 
t=0, the material is isotropic and possesses a center of symmetry, / must be 
expressible as a polynomial in the elements of an integrity basis, Los the full 
eS group of Rs eee i of the w+v(N+2) vectors Yi (« =14, 2, 
a VA (ee tee 9) are Oh re ee POP 26) 4 S20) JN andthe 
N+2 symmetric sensors G,q and Poa nabs & 4% 2, set NY sine the notation 
W_, V, and w”) for the vectors Yu" Vi and os Sd add a and the notation 
G ond yp for the symmetric Teatiiees |G, 4|| and 
elements which are non-linear in any one of the vectors W,, such an integrity 
basis, although not in general an irreducible, one, is 


VII, ( een a ie eV A 9 
OP TER (15 2,513, 0) 1, oy nee et ee 
( 4 
) 


and rit. 


where II are all possible matrix products, including the unit matrix, formed 
fromthemmatrices Gand=ve.( Pen, 4, oes). 

We denote the uw" order tensor u; ;,;,, by wu and the matrix ||éx,;/OX,|| by e. 
Also, we denote the transpose of ¢ by ec’. Then, bearing in mind the relation 
(4.9) and the definition of W, contained in (4.7), we see that u must be expressible 
as a St of outer products of order ue formed from vectors of the types V, IIc’ 
and w{”) IIe’ and the matrices eII c’, the Seen of ee terms being poly- 
nomials in invariants of the types V, I ¥,, w! PTL wy’, wl?) II V, (a, B=1, 2,...,%3 
PEO OAR? me) AG Alek 


Consequently, we must have 


Re a ler eS a 
Uj; Ve... Dios 3 (5.2) 


where ©. , 4, 18 ie ues of outer products of order «« formed from vectors of the 
types V, II and w!”) TI and matrices of the type II, the coefficients of these ake 
products being polynomials in invariants of the types V,II ¥,, wl?) TI wi? 
V, Tw!” and trII. We note from (4.10) that matrices of the type Tr and eecihe 
of the types V,II and w! II may be expressed in the forms 


oe t 
Tl fife ley ete on Tee 

00 0 

WIIG Gia Gt Glee) ene eee 
iy ye t 

VT Vf fh Pyne, hep 

00 0 3 

X ATS |G, Gale iG ta) sesey Gl te) | a ine adie 63) 

Us Y, t 

w= J ffx to, ty Tey Te) X 


x V,,(%) *[G, G (ty), G (a), ..., G(tg)] dtd ty... dtp, 
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where II* is a matrix product formed from the matrices G, G(t),..-, G(tp) 
which is linear in each of the matrices G(t,),..., G@(tp), and y(t, %, Ts, .--, Tr) 
is a continuous function of its arguments. 

In a similar manner we can express invariants of the types trII, V,ILV,, 
w)?) TI w{?) and V, Tw” in the forms 


gti t 
ta Allies Seale Tita Ast 


aS G Gia) G (2). te) | 2718 To- 12d Te, 


tt t 
VOTE, = Jal Peal Vit, ty, Uo, Ue) x 
00 0 


x VI*(G, G (1), @(%), ..., @(te)] Wdtdty...d tp, 
i et t (5.4) 
w$ Tea = J SoS 0b to. Ts Tee TRH) x 


x V,(r,) II* [G, G(T), G(T), ..., @(TR)] Vp (ta41) dd... dtRdt R44, 


aT a 
Vw)” BRERA AG T,Tg,+++, TR, TR41) X 
x V,II* [G, G(t,), G(r), ..., G (te) | (tri) day... dtpdtp +4. 


Consequently w; ;,..;, may be expressed in the form (5.2), where ® ;,__ ,,, may 
be expressed as a sum of outer products of order « formed from expressions of 
the types given on the right-hand sides of equations (5.3), with coefficients which 
are polynomials in expressions of the types given on the right-hand sides of 
equations (5.4). 

In the expressions (5.3) and (5.4), we have so far allowed II* to be any matrix 
products formed from G, G(z,),...,G(tz) which are linear in each of the 
matrices G(t,),..., G(tp). We may in fact limit the forms taken by these matrix 
products to a considerable extent and this will be done in a paper which is at 
present in course of preparation. 


6. A class of physical problems for an undeformed body 
If the body undergoes no deformation, the constitutive equation (2.2) takes 


the form ; 


Sa)]. (6.1) 


e> ( 
Uj iy...ip SG ogi los 
-— 


We shall assume that the functionals F; ;,_;,, are continuous functionals of v(t) 
over the compact aggregate of these functions which are continuous in the range 
Q<tSXt and are particularly dependent on the values of these functions at t =?. 

Now, suppose we limit ourselves to the class of problems in which each of 
the vectors v'%)(t) varies in some prescribed manner with time, so that 


P 
Hos CARO (6.2) 


where the functions /(t) are specified and the vectors wie) are independent of 
time. Then, introducing (6.2) into (6.1), we see that 


Uhl ae] iy 3 ty (w, t), (6.3) 
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where fj, ;,...;, are ordinary functions of the indicated variables, determined by 
the functionals J; ;,.;, and the functions {(z). We shall assume for simplicity 
that fj,;,...;, are polynomial functions of wi), or may be represented with suf- 
ficient accuracy by such polynomials. 

By taking 0x; (t)/@X;=G;,(t) =6;; in (3.11) and bearing in mind the defini- 
tion of V(x) given in (3.5), we see that the assumption that the constitutive 
equation is unaltered by a simultaneous rotation of the reference system and 
the physical system imposes no restriction on the constitutive equation in the 
case when the body does not undergo deformation*. Consequently, we may pass 
immediately to the consideration of the effect of symmetry of the material in 
restricting the form of the constitutive equation. 

These restrictions arise, as previously discussed, from the fact that the con- 
stitutive equation must be form-invariant under the group of transformations 
{st relating a complete set of equivalent coordinate systems in the material. 
Then, following an argument similar to that employed in § 4 and in previous 
papers, we obtain ok EF 

Ui is...ip By Oy... Oy’ (6.4) 
vat % Le 


where F is a polynomial scalar invariant under the group of transformations {s} 
of the y vectors w® (dy 2555 1e) anc the ss vectors pi?) (8.219 2 oc sh 
multilinear in the vectors y!*) and is also a single-valued function of ¢. Apart 
from the dependence on ?, this is the form which would be obtained for 1, j,i, 
if we had considered a constitutive equation relating a tensor mj, ;, __;,, and » vec- 
tors wie) for a material of the same symmetry, for which the relation is not time- 
dependent. Such constitutive equations arise as special cases, in which the 
deformation is taken as zero, of the constitutive equations discussed in the 
previous paper by PIPpkKIN & RIvLIN [6]. In that paper, the restrictions on the 
form of the constitutive equations imposed by the symmetry of the material are 
obtained in explicit form for a number of types of symmetry. These results may 
be readily applied to the case under consideration at present by introducing the 
fact that the body is undeformed and inserting dependence of UE ei ODL 
For example, in the case when « =2, and the material considered is isotropic 
and possesses a center of symmetry, it is seen from the results obtained in [6] 
that if we take as our starting point a constitutive equation of the form 


55 = fej (W”), (6.5) 


where /;; are polynomials, then «;; must be expressible in the form 


v 
u;,= 2, Ape w+ B Oi; (6.6) 
a, B=1 

* If we examine the effect of this assumption on the form of equation (6.3) 
directly, we find that it is equivalent to the restriction that the material is initially 
isotropic, but does not necessarily possess a center of symmetry. This apparent 
paradox vanishes if we realise that when the body and associated vector fields are 
rotated, the constitutive equation in the unrotated coordinate system » will not have 
the same form as it does originally, unless we include, as dependent variables in the 
initial equation, three mutually perpendicular vectors which rotate with the body 
and serve to define the directions of the symmetry axes relative to the reference 
system. If the deformation gradients occur in the constitutive equation, they provide 

these vectors. i 
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where A,, and B are polynomials in wi) wl?) It follows immediately that if 
we take (6.3) as our starting point, then u;; is expressible in the form (6.6), 
where A,, are polynomials in wi) wl?) and functions of ¢. It may be noted that 
the same result is obtainable from the results of §5 of the present paper by 
taking G(r) =7 and introducing the expressions (6.2) for vi”. 


7. A class of physical problems for a body which undergoes 
small deformations 


If we take as our starting point the constitutive equation (2.2), then we 
have seen in §3 that the assumption that it is unaltered by a simultaneous 
rotation of the reference system and the physical system implies that it can be 
expressed in the form (3.11). 

We shall now make the assumption that the deformation to which the body 
is subjected is sufficiently small, so that the displacement gradients may be 
neglected in comparison with unity. Let u;(t) be the displacement components 
in the coordinate system x undergone at time t by a generic particle of the body 
which is at X, prior to time t =O. Then, 


4;(t) = X;+ u,(z). (7.1) 


Introducing (7.1) into the definition (3.5) of G,,(7) and neglecting terms of second 
degree in the displacement gradients in comparison with those of the first degree, 
we have 


Gp q(t) — Opq= 2ep,(T), (7.2) 
where 
OUy (7) du, (T) 
toa) = (ax, tax, (75) 


Also, from (3.5), we see that neglecting du,(t)/0X, in comparison with unity, 


we have 
Vea) a). (7.4) 


We now limit ourselves to the class of problems in which the displacement 
vector u,(t) and each of the vectors v(t) varies in some prescribed manner 
with time so that 
Uy (T) (t) U, 


75 
f(r) = 70%) wf Cie 


and 


, 


where the functions g(t) and /™(z) are specified and the vectors U, and wl) are 
independent of time. Then, introducing the first of the relations (7.5) into (7.3) 
and the relation so obtained into (7.2), we have 


Gpq(t) — bpq= 28 (t) Eng, (7.6) 
where 
51 (ou au, ry 
Eng elias ! tt (7.7) 


Using (7.4), the second of the relations (7.5) and the relation (7.6) in the func- 
tional ®, ;,_;, occurring in (3.11), we see that it can be expressed as an ordinary 


ra ie ; z 


nction of w) and E,,, and of é, thus: | 


t 


Dis. in|Go ql)» VoD) | = firts.--iu 


Re t=0 

bis % ‘ j A 5 
i Introducing this result into (3.11), using (7.1) and neglecting terms of the first ~~ 
- degree in 0u,/0X, in comparison with unity, we obtain | . 
a ae te He in ale TOA ts (7-9) 
a in 


iS 


_ The restrictions imposed on the form of this constitutive equation by sym- 
metry of the material may be developed in a manner similar to that employed — 
‘in the previous cases discussed. If the group of transformations characterising 
the symmetry of the material is {s}, then 
Uist = i 
dete By OW yl’ 


(7.10) 


_ where F is a polynomial scalar invariant under the group of transformations {s} 
: ye. ‘of the » vectors w (= 15 22 ioe, yy Cheat eCOrs yi?) (B=1, 2.8.4, 8) and the 
Y symmetric second-order tensor E,,; F is multilinear in the vectors wy?) (B=41,. 
..., 4) and is also a single-valued function of ¢. 
For particular symmetries, we can again obtain explicit results, from those 
a [6] in the case when no time dependence is involved, 
wy _ by making the approximation that 0u,/0X, is small compared with unity and 
_ introducing the dependence on ¢. 
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